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Analysis  of  semiconductor  device  noise  in  the  past  has  generally  been  separated 
into  two  parts.  In  the  first  part  one  tries  to  understand  which  type  of  noise  is  the 
dominant  noise  source  and  derive  its  local  noise  strength  in  each  differential  volume. 
The  second  part  deals  with  deriving  the  transfer  function,  either  the  trans-conductance 
or  the  trans-admittance  (depending  on  the  circuit  bias  condition  at  the  contact),  to 
couple  the  local  noise  strength  from  each  differential  volume  to  the  external  contact. 
Integrating  the  product  of  the  local  noise  strength  and  the  transfer  function  over  the 
device  volume  will  produce  the  overall  noise  at  the  specified  terminal.  This  type  of 
analysis  can  only  be  done  to  a certain  extent,  since  deriving  the  analytical  expression 
of  the  transfer  function  is  limited  to  only  a few  simple  cases  and  device  structures. 

To  expand  the  capability  of  predicting  device  noise  through  this  type  of  analysis, 
we  developed  a partial  differential  equation  (PDE)  -based,  generalized-scheme  two- 
dimensional  semiconductor  device-simulation  tool  to  perform  the  transfer  function 
calculation. 

The  purpose  of  this  research  is  to  use  a computer-aided  device  simulation  tool 
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to  analyze  and  simulate  low  frequency  semiconductor  device  noise.  To  investigate 
the  low-frequency  noise,  the  bulk  generation-recombination  noise  mechanism  was 
implemented.  Since  this  type  of  noise  can  only  explain  Lorentzian  spectra  observed  in 
resistors  but  not  the  l//-like  noise  in  MOS  devices,  interface  generation-recombination 
and  the  McWhorter-type  oxide  trapping  noise  mechanisms  were  added.  The  simulation 
results  verified  that  oxide  trapping  noise  (caused  by  the  carrier  fluctuations  in  the 
channel  because  of  carrier  tunneling  between  the  trap  centers  in  the  oxide  and  at  the 
interface)  is  indeed  the  source  of  the  1 If  noise  observed  in  MOS  devices. 

One  application  of  this  tool  is  noise  defect  spectroscopy  by  modeling  the  defect 
density  distribution  in  the  oxide  through  inverse  engineering.  By  adjusting  the  defect 
density  in  the  bulk  and  in  the  oxide  to  match  the  noise  simulations  with  their  corre- 
sponding measured  data,  one  can  obtain  the  defect  density  profile  of  the  device.  In 
addition,  this  tool  can  be  used  to  predict  the  maximum  amount  of  defects  allowed  to 
guarantee  excess  noise-free  device  operation. 


Xll 


CHAPTER  1 
INTRODUCTION 

Noise  performance  of  today’s  semiconductor  devices  has  become  an  important 
factor  in  the  design  of  analog  and  VLSI  circuits.  As  the  devices  continue  to  shrink 
in  size,  voltage  and  current  signals  are  scaled  down  to  accommodate  the  changes. 
Because  the  noise  floor  does  not  scale  proportionally,  the  signal-to-noise  ratio  reduces 
as  a result.  A reduction  of  noise  can  only  be  achieved  with  better  understanding  of 
noise  mechanisms  and  noise  contributions  from  each  physical  region  of  the  device. 

Once  the  critical  regions  of  noise  are  identified,  proper  process  steps  can  be  applied 
during  the  manufacturing  processes  to  lower  the  noise  sources  in  those  regions,  in  order 
to  obtain  the  reduction  of  overall  noise  desired  at  the  device  terminals. 

1.1  Noise  Sources  in  Semiconductor  Devices 

The  microscopic  noise  mechanisms  within  semiconductor  devices  can  be  classified 
into  two  categories.  They  are  carrier  velocity  fluctuations  and  number  fluctuations. 
Carrier  velocity  fluctuations  are  generally  caused  by  the  random  scattering  of  carriers 
in  the  devices.  Their  magnitude  is  determined  by  the  carrier  diffusivity.  They  appear 
as  a white  noise,  for  which  the  noise  magnitude  is  constant  as  a function  of  frequency, 
when  observed  in  a spectral  window,  as  shown  in  Figure  1-1.  This  mechanism 
becomes  dominant  when  noise  is  measured  at  high  frequencies.  It  also  sets  the  lower 
limit  of  the  noise  at  low  frequencies.  Examples  of  this  type  of  noise  mechanism  are 
the  thermal/Nyquist  noise  from  resistive  components,  the  shot  noise  in  device  junctions, 
and  the  diffusion  noise  from  the  channels  in  MOS  structures. 

Carrier  number  fluctuations  are  caused  by  random  carrier  transitions  between 
continuum  states  and  localized  defect  states.  Fluctuations  in  the  number  of  electrons 
and  holes  in  the  conduction  and  valence  bands,  respectively,  are  attributed  to  the 
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trapping-detrapping  of  carriers  to  either  bulk  or  surface  defect  centers.  Each  trap 
center  is  characterized  by  its  energy  level,  trap  density,  and  carrier  capture  coefficients. 
The  trapping  noise  appears  to  be  constant  at  low  frequency,  and  then  rolls  off  with 
a slope  of  1 If2  when  it  passes  its  comer  frequency,  which  is  determined  by  the  trap 
characteristic  time,  as  shown  in  Figure  1-1.  This  phenomenon  can  be  seen  in  resistor 
and  JFET  structures. 

Low  frequency  1 //  noise  is  typically  found  in  all  active  devices,  and  is  also  shown 
in  Figure  1-1.  There  are  debates  on  the  causes  of  such  phenomenon.  One  proposed 
by  Hooge  [1]  suggests  that  it  is  from  carrier  diffusivity/mobility  fluctuations  due 
to  the  electron-lattice  interaction,  such  as  acoustic  scattering  events.  Furthermore, 
the  noise  magnitude  is  shown  not  to  scale  significantly  with  impurity  concentration. 
Measurement  data  obtained  by  the  author  [2]  on  several  n+-p  gated  diodes  showed 
otherwise. 

The  other  possible  explanation  of  1//  noise  was  first  proposed  by  McWhorter  [3], 
and  assumes  that  superposition  of  a large  number  of  Lorentzian  spectra  produced  by 
generation-recombination  (g-r)  mechanisms  would  generate  a l//-like  noise  spectral 
density  at  low  frequencies,  as  shown  in  Figure  1 -2.  Such  a process  can  occur  at  the 
trap  centers  in  the  bulk  and  at  the  surface.  Analytical  derivations  of  this  McWhorter- 
type  g-r  noise  in  the  last  half  century  all  pointed  out  that  this  model  applies  well  to 
the  trapping  and  detrapping  of  carriers  at  the  trap  centers  in  the  oxide,  since  only  this 
process  would  have  a characteristic  time  that  is  large  enough  to  distribute  Lorentzian 
components  at  low  frequency. 

1.2  PDE-Based  Device  and  Noise  Simulator 

Partial  Differential  Equation  (PDE),  physics-based  device  simulators  have  quickly 
gained  popularity  in  recent  years  because  of  their  advantages  over  other  types  of 
device  simulators  and  design  tools.  They  can  provide  fast  and  accurate  results.  For 
example,  they  require  fewer  computations  than  Monte  Carlo  simulators.  In  addition, 
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PDE-based  device  simulators  are  easily  adapted  to  device  model  updates.  Add-on 
software  packages  (such  as,  in  our  case,  a noise  simulator)  can  be  created  and  drawn 
from  already  produced  simulator  parameters.  Furthermore,  this  type  of  simulator  can 
be  made  to  extract  data  directly  from  a process  simulator.  This  combination  is  a cost- 
effective  tool  for  device  optimization  since  the  simulators  are  parameter  compatible. 

One  of  the  main  goals  of  this  study  was  to  implement  known  noise  mechanisms 
into  a general-purpose,  multi-dimensional,  PDE-based  device  simulation  tool,  and 
then  use  it  to  explain  different  noise  behaviors  observed  in  various  device  structures. 
This  approach  has  been  explored  by  Bonani  and  others  [4].  His  work  focused  on 
the  noise  mechanisms  mainly  observed  in  the  bulk  region  of  semiconductor  devices, 
which  are  velocity  fluctuation  noise  and  direct  band-to-band  generation-recombination 
(g-r)  noise.  He  later  extended  his  focus  to  cover  g-r  noise  from  carrier  transitions 
between  the  continuum  states  and  the  defect  trap  centers  in  the  forbidden  energy 
band-gap  [5],  around  the  same  time  we  implemented  the  same  mechanism  to  our 
noise  simulation  tool  [6].  More  on  this  topic  is  presented  in  the  next  chapter.  As 
expected  from  the  analytical  theory,  velocity  fluctuation  noise  simulation  results  in 
a white  noise  spectrum  in  the  spectral  domain,  and  bulk  g-r  noise  simulation  gives 
a Lorentzian  shape  spectrum  for  resistor-like  device  structures.  Neither  mechanism 
could  explain  the  l//-like  noise  observed  in  bipolar  junction  transistors  (BJT)  or  metal 
oxide  semiconductor  field-effect  transistors  (MOSFET).  As  the  topic  of  this  research 
states,  our  objective  is  to  use  noise  simulation  to  analyze  low-frequency  noise  in 
semiconductor  devices,  including  l//-like  noise.  In  order  to  achieve  this  objective, 
we  choose  FLOODS  (FLorida  Object-Oriented  Device  Simulator),  a PDE-based 
generalized  box  scheme  device  simulator  capable  of  doing  two-dimensional  dc,  ac,  and 
transient  simulations,  as  the  tool.  It  was  first  expanded  with  velocity  fluctuation  [7] 
and  bulk  g-r  noise  post-processor  noise  capability  [6],  as  Bonani  did  with  the  device 
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simulator  of  his  choice,  then  surface  noise  mechanisms  were  added  to  account  for  the 
l//-like  noise  observed  at  low  frequencies. 

1.2.1  Theory  of  Noise  Simulation 

Noise  implementation  is  based  on  the  impedance  field  method  (IFM)  [8].  Two 
elements  are  required  to  calculate  noise  at  the  contact  terminals.  One  is  the  magnitude 
of  the  microscopic  noise  sources  at  each  grid  point,  represented  by  the  added  Langevin 
noise  terms  in  the  PDE  equations.  Depending  on  their  physics  and  carrier  interaction 
mechanisms,  some  noise  sources,  such  as  diffusion  and  bulk  g-r  noise  sources,  can  be 
calculated  using  the  steady-state  solution  of  the  available  variables  (yr,  n,  p,  J,  yU,  etc.) 
at  each  node.  On  the  other  hand,  there  are  mechanisms,  carrier  tunneling  for  example, 
that  require  variables  from  nodes  that  are  far  apart  to  describe  their  dynamics  and  noise 
properties. 

The  other  element  to  calculate  noise  is  the  transfer  function  that  couples  the 
perturbation  of  the  PDEs  from  each  grid  point  to  the  device  terminals.  This  transfer 
function  is  generally  called  the  Green’s  function.  The  Green’s  function  for  each  carrier 
type  at  each  node  is  obtained  by  perturbing  the  PDE  associated  with  the  carrier  at 
the  node,  and  then  observing  the  changes  at  the  contact  terminals.  Assume  M is  the 
number  of  PDEs  in  the  system,  and  N is  the  number  of  nodes.  To  get  the  Green’s 
functions  for  all  the  variables  and  all  the  grid  points  for  each  contact  terminal,  it  would 
take  (M  x N ) matrix  inversions.  Bonani  [4]  simplifies  this  process  by  using  Branin’s 
method  [9],  and  then  it  takes  only  one  system  matrix  inversion  to  calculate  the  Green’s 
function  for  all  the  nodes  and  PDEs  in  the  system  for  a particular  contact. 

By  implementing  the  noise  terms  to  the  electron  and  hole  continuity  equations  as 
an  example,  we  have  the  following  PDEs 

dn  1 . , ,,  . 

~77  ~ ' CAi  + <5«)  — Un  + Yn 

dt  q 


(U) 
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dt  — ^ + Up  + Yp  (1-2) 

where  £,n  £p  and  yn  ,yp  are  the  Langevin  noise  terms  for  velocity  fluctuation  and 
number  fluctuation  noise,  respectively.  For  the  number  fluctuation  noise  sources  yn 
and  yp  , the  effect  is  transferred  to  the  specific  contact  terminal  by  the  scalar  Green’s 
function  G . For  the  velocity  fluctuation  noise  t,n  and  %p  , the  effect  is  through  the 
vector  Green’s  function  G , which  is  the  gradient  of  the  scalar  Green’s  function,  since 
it  is  the  divergence  of  the  noise  source  that  is  added  to  the  dn/ dt  , dp / dt  terms.  The 
total  voltage  noise  spectral  density  at  the  contact  terminals  is  then  calculated  from 
integration  over  all  noise  sources  [10] 

S(ffl)=  £ / GaKr^dr+  £ / Gak^&dr  (1.3) 

“>J a,P=n,pQ 

where  K and  K represent  the  microscopic  noise  sources  for  number  and  velocity 
fluctuation  noise,  respectively.  By  using  the  scalar  Green’s  function  already  calculated 
by  FLOODS  as  an  example,  the  g-r  noise  spectral  voltage  density  contact  contribution 
from  each  local  g-r  noise  source  can  be  calculated,  as  shown  in  Figure  1-3. 

The  microscopic  noise  sources  K for  number  fluctuations  are  derived  from  the 
Poisson  process  and  are  characterized  as  twice  the  carrier  transition  rates.  Examples 
are  shown  in  subsequent  chapters.  The  magnitude  of  sources  K for  velocity  fluctuation 
noise,  or  so  called  diffusion  noise,  are  given  by 


K> 


-=  4 qlnDn 


(1.4) 


K> 


4q  pDp 


K, 


^n&p 


= Kp  £ -0 


(1.5) 


(1.6) 
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The  velocity-fluctuation  noise  simulation  results  are  shown  below. 

1.2.2  Velocity-Fluctuation  Noise  Examples 

Velocity-fluctuation  noise  simulations  were  carried  out  with  FLOODS  on  different 
semiconductor  device  structures.  The  following  examples  only  consider  the  three 
basic  PDEs  for  the  drift-diffusion  model:  the  Poisson,  electron,  and  hole  continuity 
equations.  Langevin  noise  sources  for  the  diffusion  noise  sources  are  specified  for  both 
the  carrier  continuity  equations.  Noise  simulation  results  are  compared  with  expected 
analytical  results.  As  stated  previously,  velocity  fluctuation  noise  is  calculated  from 
the  vector  Green’s  function,  which  is  the  gradient  of  the  scalar  Green’s  function.  Good 
agreement  between  the  simulation  and  the  theoretical  expression  is  crucial,  since  it 
validates  both  the  vector  and  scalar  Green’s  functions  generated  by  the  simulator. 

Before  proceeding  to  any  simulation  of  a more  complicated  structure,  simulations 
on  a uniformly-doped  silicon  bar  were  first  carried  out.  The  simulated  results  are  in 
good  agreement  with  the  calculated  Nyquist  noise 

Si4iff  = 4kBTRe{Y(f)}  (1.7) 

where  Y (/)  is  the  small-signal  admittance  calculated  from  small-signal  ac  simulation 
result,  kg  is  Boltzmann  constant,  and  T is  the  temperature. 

pn  Diode.  Noise  simulation  is  carried  out  on  a 2-D  diode  structure.  A cross- 
section  of  its  doping  profile  between  the  anode  and  cathode  is  shown  in  Figure  1-4.  In 
the  0.5  fim  by  0.5  fim  p-doped  emitter  region,  the  acceptor  doping  concentration  is 
set  to  5-1017  cm-3.  The  3 jum  by  3 fim  n-doped  bulk  substrate  is  specified  to  have  a 
donor  doping  concentration  of  1015  cm-3.  In  a diode  structure,  we  expect  the  velocity 
fluctuation  noise  under  forward  bias  and  low-injection  to  show  up  as  2 ql,  the  shot 
noise,  and  then  expect  it  to  deviate  from  2ql  at  high  bias  voltage  due  to  high  injection, 
where  I is  the  diode  current.  The  expected  shot  noise  2ql  is  plotted  against  the 
simulated  results,  and  is  found  to  be  in  good  agreement  under  low  injection  conditions, 
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as  shown  in  Figure  1-5.  To  investigate  the  overall  noise  contribution,  the  distributed 
current  noise  spectral  densities  under  low  and  high  injection  conditions  are  plotted  in 
Figure  1-6  and  Figure  1-7,  respectively.  As  shown  in  Figure  1-6,  in  low  injection,  noise 
produced  by  the  excess  minority  carriers  dominates.  As  the  voltage  across  the  diode 
increases,  noise  produced  by  the  majority  carriers  also  increases.  Noise  of  the  diode 
operating  in  high  injection  is  shown  in  Figure  1-7.  In  both  cases,  noise  produced  in  the 
space  charge  region  is  not  the  main  factor.  According  to  Bonani  [4],  this  phenomenon 
is  explained  by  the  changes  in  the  vector  Green’s  function.  In  low  injection,  the  vector 
Green’s  functions  for  the  majority  carriers  are  negligible,  thus  the  contributions  by  the 
majority  carriers  are  minimal.  As  the  device  operation  moves  into  high  injection  which 
causes  the  vector  Green’s  functions  to  increase  drastically,  the  noise  contribution  of 
the  majority  carriers  starts  to  overshadow  the  minority  contributions.  This  effect  was 
confirmed  by  our  simulation  results. 

npn  BJT.  The  doping  profile  of  a 2-D  bjt  structure  is  plotted  in  Figure  1-8. 

The  n-type  emitter  region,  p-type  base  region,  and  n-type  collector  substrate  are 
uniformly-doped  with  doping  concentrations  of  1020  cm-3,  2-1017  cm-3,  and  1015 
cm-3,  respectively.  The  emitter  region  is  about  0.06  jUm  deep,  and  the  base-collector 
junction  is  0.5  jUm  below  the  surface.  For  this  device  structure,  we  expect  the  noise 
at  the  base  and  collector  contacts  to  follow  the  shot  noise  expression,  like  the  results 
we  have  seen  in  diodes.  Thus  we  compare  simulated  current  noise  spectral  densities 
at  the  base  and  collector  contacts  with  2qlg  and  2 qlc,  respectively.  The  results  for 
low  injection  conditions  are  shown  in  Figure  1-9.  There  is  good  agreement  between 
the  expected  and  the  simulated  results.  The  distributed  current  noise  spectral  density 
calculated  at  the  base  contact  is  plotted  for  electrons  and  holes,  in  Figure  l-10a  and 
1 -10b,  respectively.  Similar  to  what  we  obtained  in  the  diode  case,  noise  contributions 
are  dominated  by  the  minority  carriers.  The  distributed  noise  spectrum,  shown  in 
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Figures  1-1  la  and  1-1  lb,  shows  that  the  majority  contributions  of  the  noise  at  the 
collector  contact  are  from  the  minority  carriers  at  the  base  region. 

n-channel  MOSFET.  Noise  simulation  is  carried  out  on  a 2-D  long  channel 
MOSFET.  Its  doping  profile  is  displayed  in  Figure  1-12.  Both  the  n-type  source 
and  drain  regions  have  a 0.25  fim  deep,  1018  cm-3  donor  dopant  implantation.  The 
uniformly  doped  p-type  substrate  has  a 1017  cm-3  doping  concentration.  The  diffusion 
current  noise  is  expected  to  be  4 kgTyg where  g^o  is  the  zero  drain  bias  channel 
conductance,  and  for  long  channel  devices  y varies  from  1 to  2/3  as  the  device 
operates  from  the  triode  to  saturation  mode  in  the  absence  of  hot  carrier  effects.  The 
diffusion  noise  simulation  results  are  plotted  in  Figure  1-13.  As  expected  the  diffusion 
noise  is  lowered  to  2/3  of  the  value  recorded  at  low  Vjs,  as  the  drain  bias  increases 
beyond  Vgs  — Vt.  With  Vt  = 0.932  V,  noise  spectral  density  contributions  at  the  drain 
contact  for  electrons  and  holes  are  plotted  in  Figures  l-14a  and  1-1 4b,  respectively.  In 
the  triode  mode  of  operation,  the  noise  due  to  electrons  dominates.  The  channel  is  the 
critical  region  for  producing  the  overall  diffusion  noise  observed  at  the  drain  terminal, 
as  the  simulation  shows  in  Figure  l-14a. 

1.3  Organization  and  Overview 

This  dissertation  is  divided  into  five  chapters,  including  this  introduction.  Chapter 
2 provides  the  theory  behind  the  bulk  g-r  noise.  Two  different  implementation  methods 
are  described  in  detail.  The  “direct”  method  places  the  added  trap  electron  continuity 
equation  directly  into  the  system.  Noise  sources  in  each  equation  are  coupled  out  to 
the  external  contact  through  its  Green’s  function.  The  “indirect”  method  implements 
the  effect  of  the  added  trap  level  through  back  substitution.  Thus  the  noise  sources 
in  the  added  trap  continuity  equation  are  back-substituted  to  the  Poisson,  electron, 
and  hole  continuity  equations  before  they  are  coupled  out  to  the  external  contact 
through  the  modified  scalar  Green’s  functions.  Both  methods  are  compared,  and  their 
simulation  results  are  presented  at  the  end  of  the  chapter. 
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From  the  simulation  results  shown  in  Chapter  2,  the  bulk  g-r  noise  in  active 
CMOS  and  bipolar  devices  is  eliminated  as  the  possible  cause  of  the  1//  noise 
observed  at  low  frequencies.  The  focus  is  then  moved  to  the  Si/Si02  interface  region. 
Although  the  theory  behind  the  g-r  noise  at  the  interface  is  very  similar  to  the  theory 
in  the  bulk,  the  differences  between  these  two  cases  must  be  considered.  Chapter  3 
discusses  how  the  distributed  trap  energy  at  the  interface  is  taken  into  consideration 
by  the  simulator.  In  addition,  the  conversion  from  a 3D  bulk  g-r  simulation  to  a 2D 
interface  g-r  simulation  is  also  described  in  detail.  The  chapter  ends  with  the  interface 
noise  simulation  results. 

The  simulation  results  in  Chapter  3 also  point  out  that  the  interface  trap  noise  is 
not  the  source  of  low-frequency  1 If  noise.  Thus  the  focus  of  the  research  is  pointed 
to  the  noise  produced  by  the  defects  in  the  oxide.  Chapter  4 provides  the  theory 
behind  the  oxide  trapping  noise,  which  is  caused  by  the  trapping  and  detrapping  of 
carriers  from  the  channel  to  the  traps  in  the  oxide  via  trap  centers  at  the  interface. 

The  “Capture  and  Tunneling”  model  describes  how  carriers  first  get  trapped  by  the 
fast  defect  centers  at  the  interface,  and  then  tunnel  to  the  trap  centers  in  the  oxide 
at  equal  energy  level.  Depending  on  the  spatial  position  of  the  oxide  trap  from  the 
interface,  the  characteristic  time  that  determines  the  length  of  a carrier  staying  in  the 
trap  centers  is  distributed.  Distributions  of  different  characteristic  times  from  different 
trap  locations  can  then  be  transformed  to  the  frequency  domain.  Thus  the  Lorentzian 
spectra  from  these  sources  have  their  comer  frequencies  evenly  distributed,  in  the  case 
of  a homogeneous  trap  distribution,  to  produce  l//-like  noise  at  low  frequencies.  Two 
critical  factors  that  govern  the  outcome  of  the  oxide  trapping  noise  are  the  oxide  defect 
density  distribution  and  the  carrier  tunneling  rate.  Both  factors  affect  the  magnitude 
and  the  slope  of  the  1 If  simulated  noise  results.  The  modeling  of  the  defect  density 
distribution  in  energy  and  in  the  spatial  domain  is  discussed  in  detail,  so  is  the  shape  of 
energy  barrier  for  tunneling  carriers. 


Generation-recombination  noise  associated  with  bulk  defect  levels  is  modeled 
in  a partial  differential  equation-based  device  simulator  using  the  impedance  field 
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method  to  determine  the  maximum  allowable  defect  density  that  guaranties  generation- 
recombination  noise-free  operation  at  all  bias  conditions.  In  Chapter  5,  a detailed 
comparison  of  simulated  and  analytical  noise  data  shows  that  this  density  can  be 
estimated  accurately  even  in  the  presence  of  hot  carrier  effects  and  space  charge 
injection.  In  addition  the  noise  simulations  reveal  that  the  g-r  noise  plateau  value 
reaches  a maximum  as  a function  of  bias  and  then  decreases  at  higher  bias  due  to 
space  charge  injection  and  carrier  heating  which  reduce  the  local  strength  and  the 
transfer  function  of  the  distributed  noise  sources,  respectively. 

Chapter  6 consists  of  conclusions  and  suggestions  for  future  research. 
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Figure  1-1.  Different  types  of  noise  observed  in  a noise  experiment 
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Figure  1-2.  Low-frequency  1/f-like  noise  produced  from  a superposition  of  Lorentzian 
spectra,  caused  by  the  g-r  mechanism 
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Figure  1-3.  Local  g-r  noise  strength  to  the  external  contact  through  the  scalar  Green’s 
function 
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Figure  1 -4.  Doping  profile  of  a pn  diode.  The  emitter  region  with  a peak  acceptor 
concentration  of  5-1017  cm-3  is  implanted  in  a uniformly  doped  base 
substrate  with  a 1015  cm-3  donor  concentration. 


Noise  Current  Spectrum  (A  /Hz) 


Figure  1-5.  Current  noise  spectral  density  vs.  diode  voltage.  The  analytical  and  sim- 
ulation results  are  compared.  Note  that  the  analytical  result  is  only  valid 
under  low  injection  conditions. 
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Figure  1-6.  Distributed  current  noise  spectral  density  for  both  carriers,  under  low 
injection  conditions  (y  diode  - 0-4  V) 
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Figure  1-7.  Distributed  current  noise  spectral  density  for  both  carriers,  under  high 
injection  conditions  (V diode  = 1 -4  V) 
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Figure  1-8.  Doping  profile  of  a npn  BJT  structure 
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(a) 


(b) 

Figure  1-9.  Current  noise  spectral  density  vs.  bias  current.  The  collector  is  biased 
at  1 .0  V,  while  the  base  is  swept  from  0.3  to  0.7  V.  (a)  Noise  at  the  base 
terminal,  (b)  Noise  at  the  collector  terminal. 
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(a) 


(b) 


Figure  1-10.  Distributed  current  noise  spectral  density  at  the  base  contact  under  low 
injection  condition  (’ Vce  = 1.0  V,  V^,  = 0.5  V).  (a)  Electron  component, 
(b)  Hole  component. 
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(b) 

Figure  1-11.  Distributed  current  noise  spectral  density  at  the  collector  contact  under 
low  injection  condition  (Vce  = 1.0  V,  = 0.5  V).  (a)  Electron  compo- 
nent. (b)  Hole  component. 
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Figure  1-12.  Doping  profile  of  a n-channel  MOSFET.  The  channel  length  is  about 
2 jum. 
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Figure  1-13.  Current  noise  spectral  density  vs.  drain  bias  voltage  (Vgs  = 1.5  V) 
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(a) 


Figure  1-14.  Distributed  current  noise  spectral  density  contributions  at  the  drain  con- 
tact (Vjs  = 1.5  V,  Vgs  = 0.4  V).  (a)  Electron  component,  (b)  Hole  com- 
ponent. 


CHAPTER  2 

BULK  GENERATION-RECOMBINATION  NOISE 

2.1  Introduction 

Many  defects  that  exist  in  semiconductor  devices  are  usually  unintentional.  Ele- 
ments such  as  Carbon,  Oxygen,  and  various  metals  are  introduced  to  the  wafer  during 
process  steps,  and  then  become  carrier  trap  centers  that  affect  device  performance  and 
produce  generation-recombination  (g-r)  noise.  This  g-r  noise  is  due  to  the  fluctuation 
of  the  number  of  carriers  in  the  conduction  and  valence  bands.  Carriers  get  trapped  and 
de-trapped  from  the  defect  trap  centers  in  the  forbidden  band  gap,  causing  the  number 
of  carriers  in  the  continuum  states  to  fluctuate.  This  in  turn  causes  the  current  or  the 
voltage,  depending  on  the  condition  at  the  contact  terminal,  to  fluctuate  as  well.  This 
type  of  carrier  transition  has  generally  been  characterized  by  the  Shockley-Read-Hall 
model.  In  this  model,  each  defect  center  is  specified  by  its  defect  density,  its  carrier 
capture  coefficients,  and  its  energy  position  relative  to  either  the  conductance  or  va- 
lence band.  Due  to  the  possible  varieties  of  defects  present  in  the  device,  there  might 
be  one  dominant  defect  center,  or  multiple  defect  centers  across  the  entire  active  region. 
In  addition,  defects  are  distributed  over  space  as  they  are  more  likely  to  accumulate  at 
or  near  the  surface  or  the  metallurgical  junctions. 

As  mentioned  previously,  there  are  different  factors  that  would  affect  the  outcome 
of  the  g-r  noise.  For  example,  device-operating  temperature  governs  the  Fermi  level, 
which  sets  the  trap  occupancy  factor  that  limits  the  probability  of  carrier  transitions. 

The  carrier  transition  rate  increases  as  the  Fermi  level  approaches  the  trap  energy 
position.  When  the  Fermi  level  lines  up  with  the  trap  energy  position,  the  traps  are  half 
filled,  allowing  maximum  probability  of  the  carriers  to  get  trapped  and  de-trapped  from 
the  defect  centers.  Since  the  Fermi  level  is  a function  of  the  carrier  concentration,  the 
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doping  concentration  is  another  factor  to  determine  which  trap  states  are  critical  for 
generating  the  g-r  noise.  Another  factor  that  governs  the  g-r  noise  is  the  bias  across 
the  device.  Both  the  local  noise  strength  and  the  Green’s  function,  which  transfers  the 
fluctuation  at  each  differential  volume  of  the  device  to  a specific  contact,  change  with 
the  bias.  Thus  we  also  expect  the  magnitude  of  the  g-r  noise  to  change  with  the  bias. 
For  example,  the  simple  analytical  derivation  of  the  g-r  noise  of  a bulk  resistor  will 
show  that  the  current  noise  spectral  density  is  proportional  to  / 2 or  E2,  where  I is  the 
current  and  E is  the  field  across  the  device. 

2.2  Theory 

As  mentioned  previously,  the  g-r  process  can  best  be  described  by  the  Shockley- 
Read-Hall  model,  as  shown  in  Figure  2-1.  For  illustration  purpose,  assume  there  is  one 
defect  trap  center  in  the  bandgap,  and  the  traps  are  acceptor-like.  In  addition  to  the 
three  basic  equations 

F*  = —£T-\\P-n+ND-NX-»']=  0 (2.1) 

dn  1 

= g^'Jn  -8n  + rn-Yn  ( r,t ) = 0 (2.2) 

Fp  = ~ 8p  T fp  + Yp  (rd)  = 0 (2.3) 

a trapped  electron  continuity  equation  is  added  to  describe  the  carrier  transition  at  the 
trap  level  [11] 

„ dn,  , . 

Fn,  = +gn  ~ rn  ~ gp  + rp  - Yt  ( r,t ) = 0 (2.4) 

where  y„,  yp,  and  y are  the  Langevin  noise  terms  describing  transition  rate  fluctuations. 
The  transition  rates  can  be  expressed  as 

rn  — cnn{Np  — nt)  (2.5) 
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gn  = cnn\nt  (2.6) 

rp  = cppnt  (2.7) 

gP  = cpp\(NT-nt)  (2.8) 

where  n\  and  p\  are  the  Shockley  parameters 

n\  = NcZ\/2  (VC)  ~Ncex P (~~^C)  (2-9) 

Pi  = A^v3i/2  ~^vexp  : ^ (2.10) 


with  the  assumption  that  the  trap  level  is  at  least  a few  kT  from  the  conductance  and 
valence  bands.  In  these  equations  the  defect  species  are  specified  by  three  parameters: 
trap  energy  level  (. Ej ),  trap  density  (Nt),  and  capture  coefficients  (cn,  cp). 

Two  approaches  can  be  taken  to  solve  for  y/,  n,  p,  and  nt.  One  approach,  called 
the  “direct”  approach,  is  to  place  the  four  listed  equations  into  FLOODS  and  solve  for 
the  four  variables  simultaneously.  The  other  approach,  called  the  “indirect”  approach, 
is  to  eliminate  the  added  trap  electron  continuity  equation  by  mapping  the  trapped 
electron  density  nt  and  the  noise  term  yf  from  Equation  (2.4)  into  Equation  (2.1),  (2.2), 
and  (2.3).  Choosing  the  best  method  to  implement  depends  on  how  the  simulator 
formulates  the  PDEs  to  the  system. 

In  the  previous  versions  of  FLOODS,  each  PDE  was  hard  coded  into  the  system 
using  C++  scripts.  Since  hard-coding  extra  equations  into  the  system  is  difficult  to 
implement,  it  is  easier  to  just  modify  the  pre-existing  components  of  the  Jacobian 
matrix  for  the  three  basic  equations  of  the  drift-diffusion  model.  In  this  case,  the 
second  approach  is  preferred.  Each  trap  electron  continuity  equation  is  back  substituted 
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to  the  original  three  equations  through  modifications  of  the  matrix  components.  The 
advantage  of  this  approach  is  that  no  additional  data  storage  is  needed  and  little 
programming  overhead  is  involved,  since  the  three  basic  equations  for  the  calculation 
of  t//,  n,  p are  already  implemented  in  FLOODS.  In  addition,  since  there  are  only 
three  equations  and  three  variables  to  solve,  a faster  convergence  can  be  achieved  in 
simulation. 

The  current  version  of  FLOODS  was  rewritten  to  accept  scripting  of  PDEs.  The 
user  can  now  add  or  delete  equations  in  the  system  on  demand  by  describing  the 
equations  in  the  command  scripts.  In  this  case,  the  direct  approach  is  preferred,  since 
adding  the  extra  trap  continuity  equation  to  account  for  the  added  trap  level  can  easily 
be  specified.  The  programming  overhead  for  adding  equations  to  the  system  is  not  an 
issue  anymore.  Although  the  benefits  which  come  from  the  indirect  approach  disappear 
in  this  method,  the  noise  contributions  from  different  trap  levels  can  now  be  displayed 
and  separated  individually.  Both  methods  of  implementation  are  discussed  in  detail 
below. 


2.3  “Indirect”  Approach 

The  full  ac  noise  matrix  equation  is  defined  as  [7]  [Jc]  = [5]  where  [7]  is  the 
small-signal,  frequency  domain,  Jacobian  matrix  defined  as 
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Knowing  the  terms  g,  g,  g,  g,  g,  g,  and  g allows  n,  to  be  back- 
substituted  into  the  first  three  equations  to  get 
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where  By,  Bn,  and  Bp  are  the  noise  terms  in  the  modified  Poisson’s,  electron,  and  hole 
continuity  equations,  respectively.  These  terms  are  functions  of  yt,  yn,  and  yp,  and 
equal  to 
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where  the  trap  time  constant  is  defined  as 

T,  = [c„(n  + «i)+cp(p  + pi)]_1 


(2.16) 


(2.17) 


cn  and  cp  are  the  electron  and  hole  capture  coefficients,  and  n\  and  p\  are  the  Shockley 
parameters.  The  strength  of  the  modified  noise  sources  for  GR  become 
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plus  the  other  less  important  cross-correlated  terms.  The  variables  Xn  and  Xp  are 
defined  by 

Xn  = Y (2'21) 

Xp  = y,  = -mxp  (2.22) 

The  g-r  noise  source  strengths  defined  in  [11]  are 

Kyn,Y„  = ~Sn  +2r„  = —Kynjl  (2.23) 


KyP,yP  — 2 8p  + 2 rp  — -KypJt  (2.24) 

Kyt,yt  ~ Kyn,yn  + KyP,yP  (2.25) 

Implementation  of  the  equations  above  into  FLOODS  enables  the  program  to  calculate 
the  g-r  noise  voltage  spectral  density  at  an  external  contact  from 

Nt  rap  . 

Sv,gr=  12  / GaKya,ypG*pdr  (2.26) 

»'=  1 a,p=y,n,pJvo 1 

where  Ntrap  is  the  number  of  trap  levels  in  the  device  and  G is  the  Green’s  function 
modified  to  account  for  the  correlation  between  trap  transition  rates. 

As  mentioned  previously,  to  reduce  the  overhead  of  adding  the  trapped  electron 
continuity  equation  as  the  forth  equation  to  the  system,  the  dc  value  of  n,  can  be 
analytically  calculated  and  back  substituted  into  the  first  three  equations.  Under  dc 
condition,  gn  — rn  — gp  — rp.  The  dc  value  of  nt  can  be  obtained  from  this  equality: 

(cnn  + cpPl)NT 

nt  — — 7 ^ 7 r (Z.Z/J 

cn(n  + ni)+cp  {p  + p i) 

By  substituting  this  into  Poisson’s,  electron,  and  hole  continuity  equations,  the  changes 
in  the  device  dc  characteristics  due  to  the  effect  of  traps  are  accounted  for. 
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A similar  approach  has  also  been  explored  by  Bonani  [5]  independently,  showing 
again  that  the  Langevin  noise  terms  associated  with  trap-assisted  transitions  specified  in 
the  trap  electron  continuity  equation  are  back-substituted  to  the  Poisson,  electron,  and 
hole  continuity  equations.  Thus  the  modified  Langevin  noise  terms  exist  in  all  three 
PDEs.  The  local  noise  source  strengths  from  his  derivation  are  equivalent  to  the  ones 
in  Equations  (2. 1 8)-(2.20)  above. 

2.4  “Direct”  Approach 

As  mentioned  previously,  the  full  ac  noise  matrix  equation  is  defined  as  [J\  [Jc]  = 

[b].  Since  there  is  no  back  substitution  to  eliminate  any  equation,  the  Jacobian 
matrix  [J\  in  Equation  (2.11)  is  the  exact  matrix  implemented  in  the  system.  The 
matrix  [x\  and  [b]  implemented  in  the  system  are  also  the  same  as  the  one  specified  in 
Equation  (2.12).  The  g-r  noise  strengths  are  defined  in  Equation  (2.23)  through  (2.25). 
The  g-r  noise  voltage  spectral  density  can  then  be  calculated  from 

Ntrap  „ 

$V,gr=  / GaKyctJp  &*pdr  (2.28) 

«=1  a,P=\i/,n,p,n,Jvo1 

where  Ntrap  and  G are  the  number  of  trap  levels  and  the  scalar  Green’s  function, 
respectively. 

2.5  Simulation  Verification 

Extensive  research  on  bulk  g-r  noise  has  been  done  by  Van  Vliet  and  Van  der 
Ziel,  for  either  single  or  multiple  defect  trap  systems.  As  stated  in  their  published 
work  [11,  12],  the  derivation  of  the  noise  spectrum  starts  with  the  fluctuation  in  the 
trap  electron  continuity  equation 


dAnt 

dt 


where  the  trap  time  constant  z,  is  defined  as 


(2.29) 


-l 


T,  = [C„(«  + /Ji)+Cp  (p  + Pl)] 


(2.30) 


Taking  the  Fourier  transform  of  Equation  (2.29),  the  fluctuation  spectrum  for  the 
number  of  trapped  carriers  in  a small  differential  volume  Av  equals 
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where  (D  = 2nf.  From  the  Poisson/Gaussian  law,  the  spectrum  of  yt  is  equal  to 


(2.31) 
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(2.32) 


Equation  (2.31)  becomes 
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(2.33) 


where  the  value  at  low  frequency  plateau  is 


5 Ant,AV  {f  > 0)  — 2 ( gn  + f"n  4"  gp  + f p)  T, 


(2.34) 


To  verify  this  expression,  a 2-D  n+-n-n+  bulk  resistor  is  created.  The  n-  and  n+-region 
are  doped  with  Nj  = 1015  and  = 1018  cm-3  donor  ions,  respectively.  The  length 
of  the  n-region  L is  2 ;Um,  and  the  n+-regions  are  0.5  /tm  long  on  each  side.  Its  cross- 
section  area  A is  1 /im2.  A defect  center  Ej,  0.3  eV  below  the  conduction  band  Ec,  is 
specified  in  the  simulation.  Its  capture  coefficients,  c„  and  cp,  and  defect  density  AV 
are  set  to  10~8  cm3/s  and  1012  cm-3,  respectively.  Since  Ej  is  much  closer  to  Eq  than 
Ey,  one  can  safely  assume  gn,rn  » gp,rp  and  gn  « rn.  Equation  (17)  above  is  then 
simplified  to 


SAn,,AV  = 


4r„T2 
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(2.35) 


(l  +£02t,2)  Av  (1  + £02t,2)Av 
Constrained  by  the  law  of  space  charge  neutrality,  the  fluctuation  in  the  number  of  trap 
carriers  at  the  trap  level  is  equal  to  the  fluctuation  of  carriers  in  the  continuum  states. 
Thus, 


5a nt  — 5/V 


(2.36) 
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where  N is  the  total  number  of  carriers  in  the  conduction  band.  For  number  fluctua- 
tions such  as  the  bulk  g-r  noise  in  an  ohmic  device,  the  noise  spectral  density  relation 
is 


volume  of  the  n-region,  and  N with  LANj,  the  current  spectral  density  of  the  resistor 
now  equals 


analytical  expression  of  the  current  noise  spectral  density  of  this  structure  at  low 
frequencies  becomes 


above  is  simulated.  The  current  noise  spectral  density  from  the  simulation  is  plotted 
against  the  analytical  g-r  noise  expression  (2.38),  shown  in  Figure  2-2.  The  simulated 
noise  is  in  good  agreement  with  the  theory.  The  simulated  bulk  g-r  noise  rolls  off  at 
the  calculated  comer  frequency  fc  = 1 /2nxt  « 2 MHz,  and  has  a slope  of  1 //2  beyond 
fc.  These  results  are  indeed  able  to  prove  the  validity  of  the  bulk  g-r  noise  simulation 
package  in  FLOODS.  Simulations  results  from  sweeping  the  trap  energy  level  through 
the  upper  half  of  the  energy  band  gap  are  plotted  in  Figure  2-3.  As  expected  the  noise 
spectral  density  calculated  at  / = 1 Hz  reaches  a maximum  when  the  trap  energy  level 
is  at  or  near  the  Fermi  level.  When  the  Fermi  level  equals  the  trap  energy  level,  the 
trap  is  half-occupied  by  the  carriers.  This  condition  meets  the  highest  carrier  transition 


(2.37) 


where  V and  I are  the  dc  biases  of  the  device.  Substituting  Av  with  LA,  the  entire 


(2.38) 


Since  c„  {n  + n\)  » cp(p  + p\),  r,  can  be  simplified  as  t,  « [c„(n-l-wi)]  . The 


(2.39) 


where  n ~ Nj.  For  verification  the  calculated  value  of  this  equation  is  used  against  the 
simulation  data. 


The  g-r  noise  for  the  n+-n-n+  bulk  resistor  with  the  defect  trap  center  described 


34 

probability,  thus  in  turn  maximizes  the  g-r  noise.  As  for  this  device  with  = 1015 
cm-3,  the  Fermi  level  is  calculated  to  be  0.265  eV  below  the  conduction  band. 

The  g-r  noise  simulations  on  an  old  version  of  FLOODS  with  the  “indirect” 
approach  were  compared  with  the  simulations  on  the  current  version  of  FLOODS  with 
the  “direct”  approach.  The  simulation  results  came  out  the  same  as  expected. 

Measurements  on  a multi-trap  p+-p-p+  silicon  resistor  [13]  are  also  used  to 
validate  the  noise  simulation.  The  voltage  noise  spectral  density  and  the  characteristic 
time  Tf  of  each  g-r  trap  state  are  simulated  as  a function  of  temperature.  The  purpose 
for  varying  the  temperature  was  to  shift  the  Fermi  level,  as  shown  in  Figure  2-4.  As 
the  Fermi  level  crosses  a trap  state,  the  low  frequency  g-r  noise  plateau  value  due  to 
that  trap  will  reach  a maximum,  as  shown  in  Figure  2-5.  In  this  figure,  the  open  circles 
and  the  solid  curves  represent  the  measured  and  simulated  data,  respectively.  This 
representation  will  be  used  for  the  subsequent  figures,  unless  specified  otherwise.  The 
four  maxima  suggest  there  are  four  trap  states  that  are  electrically  active.  By  varying 
the  energy  position  of  each  trap  state  and  matching  the  noise  simulation  results  with 
the  measured  data,  the  energy  position  of  each  trap  can  be  found.  By  adjusting  the  trap 
density,  the  magnitude  of  the  simulated  noise  level  can  be  matched  to  the  measurement 
data.  At  each  temperature,  the  characteristic  time  T,  of  each  g-r  trap  was  determined. 
The  results  are  shown  in  Figure  2-6.  By  adjusting  the  values  of  the  hole  capture 
coefficients,  the  simulation  could  match  the  measured  data.  Since  the  noise  spectral 
density  of  each  g-r  trap  is  affected  by  its  energy  level  position,  hole  capture  coefficient 
cp,  and  trap  concentration  Nj,  the  FLOODS  simulations  and  procedures  of  adjusting 
these  parameters  need  to  be  done  iteratively.  The  results  are  shown  in  Table  2-1. 

G-R  noise  simulations  were  carried  out  on  the  same  NMOS  structure  that  was 
used  in  the  example  for  the  diffusion  noise  in  Chapter  1,  of  which  the  channel  length  is 
2 |tm.  In  the  first  test  with  this  structure,  three  defect  trap  centers  are  specified  in  the 
bulk,  one  located  at  £,■  (specified  to  be  Acceptor-like),  one  0.3  eV  below  Eq  (specified 
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to  be  Donor-like),  and  the  other  one  0.3  eV  above  Ey  (specified  to  be  Acceptor-like). 
Allowing  defect  densities  to  be  uniformly  distributed  and  equal  to  Nj  = 1010  cm-3, 
and  the  capture  coefficients  cn  = cp  - 10-7  cm3/s  to  be  the  same  for  all  three  traps, 
the  gate  bias  of  this  device  is  first  driven  to  1 .5  V.  For  this  device  the  threshold  voltage 
extracted  from  the  dc  simulations  is  equal  to  0.932  V.  Noise  simulation  results  at  1 
Hz  are  plotted  in  Figure  2-7  with  the  drain  to  source  voltage  sweeping  from  the 
triode  to  the  saturation  mode  of  operation.  From  these  results,  we  found  the  traps  with 
the  parameters  specified  above  cannot  produce  an  overall  g-r  noise  higher  than  the 
background  diffusion  noise.  Another  observation  is  that  the  defect  centers  at  or  near 
Ei  are  most  effective  on  producing  the  bulk  g-r  noise  in  this  MOSFET  structure.  This 
suggests  the  bulk  g-r  traps  in  the  space  charge  region  play  an  important  role,  since 
this  region  is  where  the  Fermi  level  crosses  the  energy  levels  of  those  defect  centers. 
Taking  a closer  look  of  where  the  noise  is  generated,  the  distributed  current  noise 
spectral  density  of  the  bulk  g-r  noise  for  the  device  operating  in  the  linear  mode  is 
plotted  in  Figure  2-8.  We  found  that  in  addition  to  the  noise  produced  from  the  space 
charge  region  due  to  the  reason  described  above,  the  majority  of  the  bulk  g-r  noise 
contribution  is  stemming  from  the  channel  of  the  device,  due  to  the  vast  amount  of  the 
carriers  in  that  region. 

Since  the  first  test  has  shown  that  the  bulk  g-r  noise  is  mostly  coming  from  the 
trap  centers  at  or  near  £),  in  the  second  test  only  the  trap  at  Ei  is  considered.  The 
defect  density  is  again  set  to  Nj  = 1010  cm-3.  This  time  the  capture  coefficients  are 
swept  from  10~6  to  10-8  cm3/s.  The  results  are  shown  in  Figure  2-9.  They  illustrate 
that  the  bulk  g-r  noise  increases  as  the  capture  coefficient  decreases.  In  addition,  it 
has  the  potential  to  become  higher  than  the  diffusion  noise.  To  investigate  this  further, 
noise  as  a function  of  frequency  is  simulated.  The  result  is  plotted  in  Figure  2-10. 

As  this  figure  shows,  changing  the  capture  coefficient  does  not  just  change  the  low 
frequency  noise  magnitude,  but  it  also  changes  the  comer  frequency.  Even  though  the 
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bulk  g-r  noise  from  the  trap  centers  with  the  capture  coefficient  of  10-8  cm3/s  manages 
to  become  higher  than  the  diffusion  noise,  its  effect  is  not  significant  to  produce  a g-r 
bump  on  1/f-like  noise  that  has  a comer  frequency  in  the  MHz  range.  The  final  test 
focuses  again  on  the  defect  centers  at  £,.  With  the  capture  coefficients  set  to  10-7 
cm3/s,  the  defect  density  Nj  is  swept  from  109  to  1011  cm-3.  The  result  is  plotted  in 
Figure  2-11.  As  expected,  a higher  defect  density  produces  higher  bulk  g-r  noise. 

These  simulations  have  demonstrated  that  the  g-r  bulk  mechanism  in  MOS 
structures  is  not  the  dominant  noise  source  at  low  frequencies.  With  reasonable 
parameters,  it  is  hardly  able  to  overshadow  the  diffusion  noise.  For  the  traps  that 
produce  noise  higher  than  the  diffusion  noise,  their  noise  spectrum  is  usually  not 
high  enough  to  surpass  the  l//-like  noise  that  dominates  at  low  frequencies.  Yet  the 
Lorentzian  produced  by  the  bulk  g-r  noise  might  be  one  of  the  candidates  that  produces 
the  Random  Telegraph  Signal  noise  that  has  been  seen  in  ultra  short-channel  devices. 
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Table  2-1.  Information  on  each  G-R  trap 


Energy  position 
(eV) 

Hole  capture  coeffi- 
cient 
(cm3s-1) 

Concentration 

(cm-3) 

T1 

Ev  + 0.405 

6.0- 10-9 

2.0- 1011 

T2 

Ev  + 0.300 

4.0- 10-9 

1.1-1012 

T3 

Ev  + 0.155 

9.0-10-10 

2.0- 1010 

T4 

Ev  + 0.120 

2.5-10-10 

1.0-109 
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Figure  2-1.  Carrier  transitions  between  the  continuum  states  and  localized  defect  states 


Noise  Current  Spectrum  (A  / Hz) 
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Figure  2-2. 


Comparison  between  the  simulated  and  analytical  current  noise  spectral 
density  as  a function  of  frequency.  The  bias  across  the  device  is  set  to 


0.1  V. 


Noise  Current  Spectrum  (A  /Hz) 


40 


Figure  2-3.  Bulk  g-r  noise  simulated  as  a function  of  trap  energy  level  (/  = 1 Hz) 


Energy  Level  (eV) 


Figure  2-4.  Fermi  and  g-r  trap  energy  levels.  The  energy  levels  of  g-r  traps  are  repre- 
sented by  the  dashed  lines,  and  the  position  of  the  Fermi  level  is  shown  by 
the  solid  line. 


Sv(0)/V02  (s) 


Figure  2-5.  Low  frequency  g-r  noise  plateau  Sv( 0)/Lo2  versus  1 000/7" 


Figure  2-6.  Characteristic  time  xt  of  different  traps  versus  1000 /T 


Noise  Current  Sepctrum  (A  /Hz) 
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Figure  2-8.  Distributed  current  noise  spectral  density  (/  = 1 Hz)  of  the  MOS  structure 
operating  in  the  triode  mode  of  operation  (Vgs  = 1.5  V,  = 0.3  V,  and  Vt 
= 0.932  V) 


Noise  Current  Sepctrum  (A  /Hz) 
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Figure  2-9.  Noise  simulated  with  different  capture  coefficients 


Noise  Current  Sepctrum  (A  /Hz) 
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Figure  2-10.  Noise  with  different  capture  coefficients  plotted  as  a function  of  fre- 
quency 


Noise  Current  Sepctrum  (A  /Hz) 
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Figure  2-11.  Noise  simulated  with  different  defect  densities 


CHAPTER  3 

INTERFACE  GENERATION-RECOMBINATION  NOISE 

3.1  Introduction 

The  number  of  defects  at  a device  interface  is  generally  considered  to  be  large 
when  compared  with  the  defects  in  the  bulk  or  in  the  oxide.  This  is  due  to  the 
termination  of  periodic  structures  and  lattice  mismatch.  In  some  cases,  the  defects  are 
also  due  to  the  presence  of  oxygen  and  other  impurities.  These  defects  can  become 
the  generation-recombination  (g-r)  centers  for  the  carriers  at  the  surface.  Sometimes 
they  also  become  the  stepping-stone  for  the  carriers  to  tunnel  into  the  defects  in  the 
oxide.  For  MOS  devices,  surface  g-r  noise  mechanisms  like  interface  trapping  noise 
and  oxide  trapping  noise  have  strong  effects  on  the  devices  overall  noise  performance, 
since  carriers  are  moving  toward  the  surface  to  form  a channel.  In  this  chapter,  the 
interface  g-r  noise  mechanism  will  be  discussed  in  detail. 

3.2  Theory 

Unlike  the  defect  traps  in  the  bulk  that  are  discrete  in  nature,  the  defect  traps  at 
the  interface  and  in  the  oxide  are  distributed  in  energy  over  the  bandgap,  due  to  the 
variety  of  defects  available  at  the  interface.  A few  techniques  can  be  used  to  measure 
the  trap  density  distributions  at  the  interface.  Among  them  are  Deep-Level  Transient 
Spectroscopy  (DLTS)  and  the  conductance  method  [14].  Although  the  DLTS  can  obtain 
good  results  for  the  traps  near  the  mid-gap  region,  it  is  less  suitable  for  traps  near  the 
band  edges.  In  comparison,  the  conductance  method  gives  more  accurate  and  reliable 
results  for  the  traps  distributed  over  the  bandgap.  An  example  of  measurement  results 
is  shown  in  Figure  3-1  [14].  In  the  past,  this  trap  density  distribution  was  modeled  by 
an  exponential  function 

Nt{E)=N0exp[S  \E-Et\]  (3.1) 
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where  Nq  is  the  measured  defect  density  at  Since  the  interface  traps  are  distributed 
in  energy,  decisions  have  to  be  made  on  which  energy  levels  needed  to  be  discretized. 
One  extreme  is  to  discretize  the  entire  bandgap  in  small  discrete  energy  intervals  A E 
to  account  for  all  the  trap  states  over  the  bandgap.  Although  this  implementation  gives 
the  most  accurate  results,  it  is  obvious  that  this  method  is  computationally  expensive, 
since  for  each  discrete  energy  interval  A E,  a trap  electron  continuity  equation  needs 
to  be  added  to  the  system  for  each  interface  node.  On  the  other  hand,  if  the  bandgap 
is  not  discretized  and  the  effect  of  interface  trap  states  is  implemented  through  the 
specification  of  some  analytical  representations,  the  simulation  results  will  not  be 
sufficiently  accurate,  since  too  many  approximations  would  have  to  be  made.  In 
addition,  the  interface  g-r  simulation  results  will  be  no  different  than  the  approximate 
analytical  results  that  have  been  shown  in  the  previous  published  works  in  this  area. 
Thus  this  approach  defeats  the  purpose  of  implementing  such  mechanism  to  the 
simulator. 

The  best  way  to  implement  would  be  to  discretize  only  at  the  energies  that  are 
critical  to  noise  calculation,  and  at  the  same  time  keep  the  number  of  discrete  energy 
intervals  A E small  to  avoid  any  computational  overhead.  In  addition,  the  effect  of 
non-discretized  trap  levels  needs  to  be  accounted  for.  Choosing  the  right  energy  level 
to  discretize  requires  some  knowledge  of  how  the  g-r  noise  works.  According  to  the 
classical  g-r  noise  theory,  the  power  spectral  density  due  to  carrier  density  fluctuation 
in  a differential  volume  AA  is  given  by 

SAAn,  = Y+a2x2N‘f‘(l  -ft)^  (3.2) 

where  N,  is  the  trap  defect  density,  z is  the  trapping  time  constant,  and  ft  is  the 
Fermi-Dirac  (FD)  distribution  that  determines  the  occupancy  of  traps.  It  can  be  shown 
from  Equation  (3.2)  that  the  probability  of  carriers  being  trapped  is  proportional  to 
the  function  ft  (1  — /,).  This  function,  plotted  in  Figure  3-2,  is  significant  only  at  the 
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electron  quasi-Fermi  level  Ef  at  the  interface  and  a few  kT  around  it.  The  function 
becomes  negligible  at  energies  away  from  E since  those  traps  are  either  empty  or 
filled,  depending  on  the  trap  energy  position.  In  addition,  by  integrating  Equation  (3.2) 
over  the  bandgap,  the  noise  spectral  density  in  AA  can  further  be  approximated  as  a 
delta  function  at  Ef  [15],  since 

fEc Nt(E)ft(\-ft)dE&  AkTNt  (Ef)  ft  (Ef)  [l -/,(£/)]  (3.3) 

JEV 

From  the  arguments  above,  the  need  for  discretizing  the  trap  energy  level  at  Ef  for 
noise  calculation  is  prominent.  In  addition,  the  approximation  of  taking  only  one 
discretized  trap  energy  level  is  valid  if  the  trap  defect  density  for  noise  calculation  at 
the  discretized  energy  level  is  approximated  by  4kTNt  (Ef). 

One  major  difference  between  the  implementations  of  the  bulk  and  the  interface 
g-r  mechanisms  is  that  for  the  interface,  the  trap  energy  level  is  not  constant  throughout 
the  entire  device,  but  it  tracks  with  the  quasi-Fermi  level  Ef  at  the  interface.  This 
means  a proper  dc  steady  state  solution  is  needed  at  each  interface  node  to  identify  its 
critical  trap  energy  level  before  any  interface  noise  simulation  can  be  taken.  Procedures 
to  perform  dc  and  noise  simulations  will  be  discussed  in  detail  below. 

3.3  DC  and  Noise  Simulations 

Although  the  traps  at  the  energy  levels  other  than  Ef  are  not  discretized,  the 
trapped  carriers  at  those  levels  induce  a surface  charge  Qjtf  at  the  interface  that  needs 
to  be  accounted  for  in  the  steady  state  condition.  If  the  effect  of  the  charge  Qnf 
is  not  included  in  the  system,  the  solutions  of  yr,  n,  and  p at  the  interface  will  not 
be  calculated  correctly.  This  would  introduce  error  to  the  value  of  Ef,  forcing  the 
simulator  to  discretize  at  a wrong  trap  energy  level  and  result  in  faulty  noise  simulation 
results.  Before  any  simulation,  the  traps  at  the  interface  first  need  to  be  characterized 
as  either  donor-like  or  acceptor-like.  Assuming  the  interface  traps  to  be  acceptor-like, 
as  shown  in  Figure  3-3,  the  interface  charge  can  then  be  calculated  as  the  product  of 
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the  electron  charge  q and  the  trap  carrier  density  nt(E)  integrated  over  the  bandgap 


where  Nt(E)  is  the  trap  density  described  in  Equation  (3.1),  and  ft(E,Ef)  is  the  FD 
distribution  function.  The  expression  Qpf  in  this  equation  can  be  re-written  as  a 
function  of  n and  p,  which  are  the  variables  that  the  system  is  solving  for.  This  is 
convenient  since  no  new  variable  is  needed.  At  the  interface,  the  carrier  transition  rates 
described  by  the  processes  a,  b,  c,  and  d in  Figure  3-2  are  also  needed  to  be  set  as  a 
function  of  the  quasi-Fermi  level  Ef,  or  a function  of  n and  p,  since  the  discretized  trap 
energy  level  E,  moves  with  Ef.  Because  E,  lies  at  Ef,  the  Shockley  densities  are  equal 
to 


The  rate  equations  characterizing  the  processes  a,  b,  c,  and  d in  the  figure  can  then  be 
written  as 


n\  = «,exp  \{Ef  — Ei)  /kT]  = n 


(3.5) 


p\  = rn exp  [(Et  - Ef)  /kT]  =p 


(3.6) 


rn,itf  = cnn(N't -nt)  =cnn  4kTN0eS(Ef  E‘) -nt 
= cnn  4 kTN0{n/rii))^kT  — nt 


(3.7) 


8n,itf  = cnn\nt  = cnnnt 


(3.8) 


rp,itf  — CpPnt 


(3.9) 


gP,itf  = Cpp\  (Nj-n,)  = CpP\  4 kTN0e^Ef  E‘) -n, 
= cpp  4 kTN0(n/rii))^kT  -nt 


(3.10) 


respectively.  Notice  that  the  defect  density  N't  at  the  discretized  energy  level  is 
approximated  as  4kTN,(Ef).  Equations  (3.5)  through  (3.10)  are  only  valid  for  the  dc 
simulations.  The  pertinent  equations  at  the  interface  are 
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(3.12) 


(3.13) 


where  yn,itf,  Yp,itf,  and  yt,itf  are  the  Langevin  noise  terms  describing  transition  rate 
fluctuations.  As  mentioned  earlier,  for  dc  simulation  the  values  of  r„iiry,  gn,itfi  fp,itf , 
gp,itf  and  Qitf  are  functions  of  n and  p.  Once  the  dc  solutions  are  obtained,  the 
Shockley  densities  become  constants  as  n\  — ndc  and  p\  — pdc-  The  transition  rates  are 
then  modified  to  the  following  for  noise  simulation 


rn,itf  = cnn(N'tdc-nt) 


(3.15) 


Sn,itf  — Cnn\nt  — Cnndcnt 


(3.16) 


r P,itf  = Cppn, 


(3.17) 


8P,itf  = cpP\  W,dc  ~fh)=  CpPdc  Kdc  ~ nt] 


(3.18) 
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For  both  the  dc  and  noise  calculations,  the  electron  and  hole  capture  coefficients  for  the 
discretized  trap  level  still  need  to  be  specified. 

3.4  Simulation  Results 

Interface  noise  simulations  were  carried  out  on  a 0.45  [im  nMOS  structure.  For 
the  interface  trap  parameters,  No  is  set  to  1011  cm_2eV_1  [14],  and  both  the  electron 
and  hole  capture  coefficients  are  set  to  10-7  cm3/s.  A set  of  dc  simulations  is  designed 
to  test  the  response  at  the  interface.  The  first  simulation  sets  the  drain,  source,  and 
substrate  terminals  to  ground,  and  simulates  with  sweeping  the  gate  bias  voltage. 

The  device  is  biased  to  act  as  a MOS  capacitor.  By  increasing  the  gate  bias  voltage, 
more  electrons  move  toward  the  surface.  This  shifts  the  Fermi  level  Ef  toward  the 
conduction  band.  At  the  same  time,  this  increases  the  defect  density  at  the  discretized 
defect  level.  As  expected,  a higher  trap  carrier  density  nt,itf  and  a higher  interface 
charge  <2,?/  are  shown  in  the  simulation  results,  plotted  in  Figure  3-4. 

The  second  simulation  sets  the  gate  bias  to  a value  above  V,  to  establish  a channel 
under  the  gate.  The  trap  carrier  density  at  the  interface  nodes  is  observed  as  the  bias 
Vds  increases.  The  results  are  plotted  in  Figure  3-5.  As  increases,  electrons  would 
move  away  from  the  surface  near  the  drain,  thus  moving  Ef  closer  to  £,•  and  causing 
the  trap  carrier  density  to  decrease.  In  the  next  test,  the  interface  noise  simulations 
are  carried  out  for  a few  bias  points.  The  preliminary  results  of  the  distributed  current 
noise  spectral  density  from  the  interface  nodes  under  the  gate  are  plotted  in  Figure  3-6. 
From  these  results,  higher  interface  noise  contributions  from  the  interface  nodes  near 
the  drain  can  be  observed.  Although  the  trap  carrier  density  for  those  nodes  is  lower 
than  at  the  interface  nodes  near  the  source,  higher  noise  contributions  must  come 
from  their  higher  scalar  Green’s  function.  Finally,  the  diffusion  and  interface  g-r 
noise  simulations  are  carried  out  simultaneously  on  the  testing  device  operating  from 
triode  to  saturation  mode.  The  simulation  results  are  shown  in  Figure  3-7.  For  the 
interface  defect  centers,  the  defect  density  is  modeled  by  with  an  U-shaped  energy 
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distribution  measured  from  either  C-V  or  conductance  data.  The  trap  energy  location 
of  the  discretized  energy  level  is  determined  from  the  location  of  the  Fermi  level 
at  the  interface.  Thus  the  only  unknown  variable  is  the  capture  coefficient.  For  the 
simulation  results  plotted  in  Figure  3-7,  the  capture  coefficients  for  the  interface  traps 
were  set  to  10-7cm3/s.  As  the  simulation  results  demonstrate,  noise  produced  from  the 
interface  traps  is  negligible  when  compared  to  the  diffusion  noise.  Simulations  were 
also  performed  by  varying  the  values  of  capture  coefficients  within  reasonable  physical 
range.  Those  results  again  have  all  shown  negligible  interface  noise.  This  phenomenon 
is  expected,  since  interface  traps  are  generally  considered  to  be  fast  trap  states.  Thus 
noise  produced  from  these  traps  has  a low  magnitude  at  low  frequencies  and  rolls  off  at 
high  corner  frequencies. 

The  conclusion  of  these  simulations  is  that  the  g-r  mechanism  at  the  interface 
alone  is  not  responsible  for  the  low-frequency  l//-like  noise  that  has  been  observed  in 
MOS  noise  measurement  data. 
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Figure  3-1.  Measured  trap  density  distribution  from  the  conductance  method  [14] 


Figure  3-2.  Energy  band  diagram  at  the  interface,  with  the  presence  of  a discretized 
trap  energy  level  at  Ef 


Figure  3-3.  Trapped  carriers  at  non-discretized  trap  levels  induce  a charge  Qi,f  at  the 
interface 
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Figure  3-4.  Interface  trap  carrier  density  ntjt ^(marked  by  “x”)  and  interface  charge 
Qitf/ q (marked  by  “o”)  for  a sweep  (V^  and  Vbs  are  set  to  zero) 


Trapped  Carrier  Cone,  (cm 
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Figure  3-5.  Interface  trap  carrier  density  n,^f  for  a sweep.  Vgs  is  set  to  a value 
greater  than  Vt  and  Vhs  is  set  to  zero. 
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Figure  3-6.  Distributed  interface  current  noise  spectral  density  for  different  Vds.  Vgs  is 
set  to  a value  greater  than  Vt  and  Vt>s  is  set  to  zero. 
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Figure  3-7.  Various  current  noise  spectral  densities  for  a Vgs  sweep  (V^  = 1 .0  V) 


CHAPTER  4 

OXIDE  TRAPPING  NOISE 

4.1  Introduction 

The  l//-like  noise  in  semiconductor  devices  has  been  studied  extensively  in 
the  past.  Most  evidence  has  shown  that  its  sources  are  located  at  or  near  Si/Si02 
interfaces.  This  partially  explains  the  fact  that  l//-like  noise  is  more  dominant  in 
MOS  devices  than  in  bipolar  transistors,  due  to  MOS  structures  having  larger  surface 
areas.  According  to  McWhorter  and  others  [3,  15,  16,  17,  18],  l//-like  noise  is  the 
superposition  of  different  Lorentzian  spectra  of  which  the  corner  frequencies  spread 
across  the  low-frequency  spectral  window.  The  corresponding  characteristic  time  r 
of  each  Lorentzian  spectrum,  calculated  from  its  comer  frequency,  gives  information 
on  the  tunneling  between  defect  centers  in  the  oxide  and  at  the  interface.  Even 
though  the  carrier  tunneling  mechanism  between  the  traps  in  the  oxide  and  at  the 
interface  is  the  source  of  l//-like  noise,  only  its  macroscopic  effect,  represented  by 
a characteristic  time  T,  is  used  in  analytical  derivations.  Deriving  the  l//-like  noise 
mechanism  from  carrier  tunneling  equations  and  microscopic  g-r  trapping  noise 
sources  has  not  been  attempted,  due  to  the  complexity  of  the  problem.  In  recent  years, 
device  noise  simulation  packages  have  been  implemented  successfully  into  PDE- 
based  device  simulators.  Utilizing  the  vector  and  scalar  Green’s  functions  calculated 
by  the  simulator,  the  noise  at  device  terminals  from  each  local  noise  source  can  be 
calculated.  Simulations  of  diffusion,  Hooge  1//,  and  bulk  g-r  noise  have  shown  that 
none  of  these  mechanisms  can  produce  l//-like  noise.  The  bulk  g-r  mechanisms  are 
fast  processes  with  characteristic  time  constants  too  small  to  produce  Lorentzian  spectra 
at  low  frequencies.  Surface  noise,  or  so-called  oxide  trapping  noise,  has  not  been 
implemented  until  recently.  The  difficulty  with  the  implementation  of  carrier  tunneling 
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mechanisms  and  oxide  trapping  noise  is  that  carrier  tunneling  is  not  just  dependent  on 
variables  at  local  grid  points  and  adjacent  points,  but  also  depends  on  variables  at  grid 
points  that  are  distant.  The  noise  simulator  needs  to  have  a database  to  keep  track  of 
the  link  between  each  oxide  point  and  its  corresponding  interface  point.  Recently,  the 
earner-tunneling  mechanism  and  its  associated  noise  was  successfully  implemented  in 
FLOODS  (FLorida  Object-Oriented  Device  Simulator)  [19,  20]  and  the  oxide  trapping 
noise  simulation  results  show  1/ /-like  characteristics  as  will  be  shown  below. 

4.2  Theory  and  Implementation 

As  stated  previously,  McWhorter’s  theory  has  been  used  frequently  to  explain  the 
1//  noise  in  MOS  devices.  Christensson  [15]  and  others  have  proposed  a “Tunneling 
and  Capture”  model,  in  which  (i)  the  carriers  first  tunnel  to  a distance  from  the 
continuum  states  at  the  interface  into  the  oxide,  and  then  (ii)  are  captured  by  the 
traps  that  have  the  same  energy  level  as  Ef.  This  model  is  illustrated  in  Figure  4-1. 

As  indicated  in  the  figure,  the  electron  needs  to  lose  about  0.5  eV  of  energy  to  be 
trapped  at  the  defect  level.  The  validity  of  this  model  was  later  questioned,  since 
measurement  data  did  not  support  energy  dissipation  in  the  oxide  [21].  The  “Capture 
and  Tunneling”  model  was  later  proposed  by  Fu  and  Sah  [17]  as  the  alternative  for  the 
previous  proposed  model.  In  this  model,  the  carriers  first  get  trapped  by  the  fast  defect 
centers  at  the  interface  through  the  Shockley-Read-Hall  (SRH)  process  and  then  tunnel 
into  the  traps  in  the  oxide  at  the  same  energy  level,  as  illustrated  in  Figure  4-1.  Thus 
the  earners  do  not  dissipate  energy  in  the  oxide.  The  continuous  trap  distribution  in 
energy  over  the  bandgap  at  the  interface  is  required  to  be  the  pre-existing  condition. 

It  is  validated  by  measurement  data,  since  large  varieties  of  defects  at  the  interface  is 
a common  characteristic  of  semiconductor  devices.  From  this  point  forward,  Fu  and 
Sah’s  “Capture  and  Tunneling”  model  is  used  as  the  fundamental  basis  of  the  carrier 
trapping  mechanism  in  the  oxide. 
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4.2.1  “Capture  and  Tdnneling”  Model 

Figure  4-1  shows  the  energy  band  diagram  that  illustrate  the  possible  carrier 
transition  processes  at  the  interface  and  in  the  oxide  [17].  Carriers  in  the  continuum 
states  can  get  trapped  (process  a & c)  or  de-trapped  (process  b & d)  from  the  localized 
defect  centers  at  the  interface,  through  the  SRH  process.  The  trapped  carriers  can  also 
communicate  with  defect  centers  that  have  the  same  energy  level  in  the  oxide  through 
tunneling  (process  e & f).  Since  the  trap  centers  at  the  interface  and  in  the  oxide  are 
distributed  in  energy  over  the  bandgap,  discretization  in  energy  is  needed.  In  addition, 
both  spatial  and  energy  distributions  of  interface  and  oxide  traps  are  needed.  Measured 
data  in  the  past  [22]  generally  describes  the  trap  distribution  to  be  U-shaped  in  energy. 
Thus,  it  is  modeled  as  [14] 


Nt(E)  =N0 exp  [£  \E  — £,-|]  (4.1) 

for  simplicity.  The  term  £)  is  the  intrinsic  energy  level.  Discretization  in  small  energy 
intervals  would  require  a large  number  of  trap  electron  continuity  equations  specified  to 
the  system.  Thus  carefully  choosing  the  discretized  energy  level  for  noise  simulation 
is  critical.  According  to  the  g-r  noise  theory,  the  power  spectral  density  due  to  a 
fluctuation  in  the  number  of  trapped  carriers  nt  in  a differential  volume  AA,  either  from 
the  trap  centers  at  the  interface  or  in  the  oxide,  is  [15,  18,  23] 

SAAn,  = TT7^2N‘f‘  ( Wr)  AA  (4.2) 

where  T is  the  trapping  time  constant  and  f,  is  the  Fermi-Dirac  distribution  function. 
Plotting  the  function  ft  (1  — ft),  which  is  a strong  function  of  energy  E,  will  display 
the  critical  energy  levels  where  the  contributions  of  the  g-r  noise  are  mostly  coming 
from.  As  shown  in  Figure  4-1,  this  function  is  most  significant  at  the  electron  quasi- 
Fermi  level  at  the  interface  Ef, v-.  By  discretizing  only  the  trap  level  at  Efnj,  the 
defect  density  in  Eq.  (4.2)  is  then  modeled  as  N[  = 4kTN,  (E/n,/),  so  that  the  function 
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Ntft  (1  — ft)  can  be  represented  as  a delta  function  at  Ejn  i [15].  Since  the  discretized 
trap  energy  level  Et  moves  with  Ef„i,  the  Shockley  densities  are  equal  to 

n\  = n,exp  [(£/-£,)  / kT ] = n (4.3) 

Pi  = Hi  exp  [( Ei-Ef ) /kT\=p  (4.4) 


The  rate  equations  characterizing  the  processes  a,  b,  c,  and  d in  Figure  4-1  then 
become 


rn,itf  = cnn  (N't  -nt)  = c„n  4kTNae % (Ef  E‘)  - n, 


cnn 


4kTN0  («/«,)' *kr-nt 


(4.5) 


8n,itf  — cnn\nt  — cnnnt 


(4.6) 


rp,itf  = cppnt 

gp,itf  = cpp\  ( N't  - nt)  = cpp\  4kTN0e^Ef~E‘)  - n, 


cpp 


4kTN0  {n/ntfkT  - nt 


(4.7) 

(4.8) 


respectively.  The  basic  equations  for  the  discretized  energy  level  at  the  interface 
become 


F*  = ~li?-\\P-n+Nt-Ni- &r  ( Ef ) n,Mf\  = 0 

dn  \ , , , . 

Ffl  ~ fff  ~ - ^ ■ Jn  + yn,itf  — 8n,itf)  ~ Yn.il  f \rd,Ef)  =0 

Fp  = + - V • Jp  + (rP}itf  — gpjtf ) — Yp.it  f (r,t,Ef)  = 0 


Fn,,itf  — ^ + (gn.itf  rn,itf)  ( gp.itf  rp,itf ) 


Yt,itf  ifdiEf)  (gox  T P0x)  Yox{fd  iEf)  — 0 


(4.9) 

(4.10) 

(4.11) 

(4.12) 


where  Yn.it  Yp.it f,  and  Yt,uf  are  the  Langevin  noise  terms  describing  transition  rate 
fluctuations  at  the  interface,  and  Yox  is  the  Langevin  noise  terms  associated  with 
transition  rate  fluctuations  through  carrier  tunneling.  The  terms  gox  and  rox  are  the 
tunneling  rates  illustrated  as  processes  e and  f in  Figure  4-1,  respectively.  These 
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tunneling  rates  are  functions  of  the  tunneling  coefficient  T,  the  trap  carrier  density  of 
the  defect  centers  in  the  oxide  nt)OX,  and  the  trap  carrier  density  of  the  defect  centers  at 
the  corresponding  interface  node  nt^,f 

rox  (■*/)  — Tntitf  {Nt,ox  — nt,ox ) Ax  (4.13) 

Sox  (xi)  = T nt,itf ) nt,o; cAx  (4.14) 

Notice  that  these  rate  equations  are  not  just  dependent  on  the  local  variables  at  the 
oxide  node,  but  they  are  also  dependent  on  the  variables  at  the  interface.  From  a 
simple  classical  tunneling  model  [15,  21,  14,  24]  with  a constant  barrier  height,  the 
tunneling  coefficient  is  given  by 


T = Tq  exp  (ad) 


(4.15) 


where  a = —2yj2moxtyitf/h2  is  the  attenuating  constant  of  the  wave  function  in  the 
oxide  associated  with  the  barrier  height  0Z/ f at  the  interface.  In  addition  to  the  basic 
equations  above,  the  equations  to  be  simulated  for  an  oxide  node  at  position  x,-  are  the 
Poisson  and  trap  electron  continuity  equations 


r _ 4 

Qox  ( xi  > Ef ) 

Vx‘  dr 2 e 

q 

dn0X(xi,Ef 

r°x  — , 

rox  T Sox 

- n0X(xi,Ef ) 


= 0 


Vox  (xi  dt  Ef)  — 0 


(4.16) 

(4.17) 


4.2.2  Modeling  of  the  Defect  Density  in  Oxide 

One  of  the  main  goals  of  this  research  work  is  to  characterize  the  defect  density 
distribution  in  the  oxide  region  of  semiconductor  devices.  Knowledge  of  the  defect 
density  distribution  is  critical  for  the  operation  of  a MOSFET,  since  defects  play 
important  role  in  gate  leakage  current  and  act  as  the  source  of  low  frequency  noise. 
Different  models  for  the  oxide  defect  density  can  be  applied  to  the  simulator  to 
compare  the  simulation  data  with  the  measured  noise  results.  Thus  an  optimal  defect 
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density  profile  can  be  obtained  once  a match  between  simulations  and  measurements  is 
found  via  reverse  engineering. 

Research  in  the  past  speculates  that  the  defect  density  in  the  oxide  varies  in 
position  and  in  energy,  based  on  observations  such  as  the  ones  listed  below: 

• As  stated  previously,  the  measured  interface  trap  density  is  found  to  have  a U- 
shape  in  energy.  For  the  same  reason  that  more  defect  centers  are  found  near  the 
continuum  states,  one  can  assume  that  the  traps  in  the  oxide  have  a U-shape  in 
energy  distribution  also. 

• From  a simple  classical  tunneling  model,  the  tunneling  rate  is  an  exponential 
function  of  the  distance  between  the  interface  and  the  oxide  traps,  as  stated  in 

Eq.  (4.15),  if  each  tunneling  carrier  sees  a constant  barrier  height.  It  can  be  shown 
analytically  that  the  tunneling  rate  calculated  with  this  simple  model  leads  to  the 
ideal  1 / / noise  spectrum,  with  the  condition  that  the  defect  density  in  the  oxide 
has  to  be  constant  in  energy  and  space.  For  this  reason,  many  in  the  past  model 
the  defect  density  in  the  oxide  with  a constant  value,  for  simplicity  in  analytical 
noise  derivations  and  for  the  modeling  of  1//  noise  observed  in  measurements. 

• Defects  are  accumulated  at  or  near  the  interface.  It  has  been  reported  from  noise 
measurements  that  the  logarithmic  slope  of  1/ /-like  noise  deviates  from  —1,  as 
the  frequency  increases  in  the  spectral  window.  This  type  of  noise  is  sometimes 
called  as  “pseudo  1//  noise”.  Again  by  modeling  the  tunneling  rate  as  a function 
of  distance  and  with  a constant  barrier  height,  the  measurement  data  supports  the 
theory  of  an  increasing  defect  density  distribution  at  the  regions  near  the  interface 
in  the  oxide.  Depending  on  the  location  of  the  defects,  for  those  that  have  a 
higher  defect  density  near  the  interface,  their  effect  would  only  be  shown  at  high 
frequencies,  since  they  are  very  close  to  the  interface  and  the  tunneling  rates  for 
the  carriers  communicating  with  them  are  fast. 

Based  on  the  observations  described  here,  Celik  and  Hsiang  [14]  came  up  with  an 
analytical  expression  for  the  oxide  defect  density  as  a function  of  energy  and  distance 
from  the  interface.  In  addition,  the  expression  accounts  for  band  bending  in  the  oxide. 
The  defect  density  is  modeled  as 


N,  (E,x)  = Afoexp 


UE„El)+*%=n±+r]x 


(4.18) 


where  q (Vg  - V)x/t0X  is  the  amount  of  band  bending  at  a distance  x away  from 
the  interface.  Based  on  their  noise  measurements  on  MOS  devices,  they  found 
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N0  = 4- 1017  cm“3eV_1,  ^ = 3.1  eV-1,  A = 62  eV-1,  and  rj  = 4.5 • 107  cm-1.  Ideally, 
the  band-bending  factor  A should  be  equal  to  unity,  since  it  is  more  physical  for  the 
defect  density  to  be  modeled  by  the  exact  amount  of  band  bending  at  each  location, 
without  this  multiplication  factor.  However,  the  published  value  of  A provides  a better 
match  between  their  measured  and  calculated  data. 

4.2.3  Modeling  of  the  Carrier  Ttinneling  Rates 

In  general,  the  carrier  tunneling  rate  is  represented  by 

T = T0  exp  (ad)  (4.19) 


where  the  attenuation  constant  of  the  wave  function  depends  on  the  shape  of  the 
barrier.  It  also  depends  on  the  effective  mass  of  the  carrier  in  the  oxide,  and  the 
barrier  height  at  the  interface  There  are  two  general  models  that  can  be  used 
to  characterize  carrier  tunneling.  The  simplest  one  is  the  square  barrier  tunneling 
(SBT)  model.  As  the  name  implies,  the  carrier  sees  a square  potential  barrier  before 
tunneling.  Of  course  this  is  only  true  under  flat-band  condition.  Due  to  the  simplicity 
of  the  SBT  model,  it  is  widely  used  as  an  approximation  for  calculating  the  tunneling 
characteristic  time  constant  in  analytical  derivations  of  noise  theory.  In  this  model,  the 
attenuation  constant  is  defined  as 


oc  — ^ 2m ox  (pit  f 


(4.20) 


Setting  mox  — 0.4 ma  as  obtained  from  BSIM’s  parameter  set  [25],  a is  then  equal  to 
1.303  • 108  cm-1  with  faf  = 4.05  eV.  In  the  actual  simulation  a will  not  be  a constant, 
since  it  is  a function  of  the  barrier  height  which  depends  on  the  location  of  the 
Fermi  level.  In  the  second  model  the  tunneling  carrier  sees  a trapezoidal  potential 
barrier,  a situation  found  in  MOSFETs.  To  calculate  the  tunneling  rate,  both  the  barrier 
heights  at  the  interface  </>;, f and  at  the  trap  location  in  the  oxide  (j)ox  are  required,  as 
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stated  in  the  following  expression 


a = - 


4 \J2mox  $?//  ~ <t>oI2 


(4.21) 


3 h (tyitf  Qox) 

As  (tyitf  — (f)ox)  approach  zero,  the  value  a in  Eq.  (4.21)  approaches  the  one  in 
Eq.  (4.20). 

As  mentioned  previously,  the  tunneling  rate  equations  specified  as  Eq.  (4.13)  and 
(4.14)  rely  on  non-local  variables.  The  effect  of  tunneling  from  any  oxide  node  to  each 
interface  node,  and  vice  versa,  needs  to  be  examined.  One  specific  behavior  of  carrier 
tunneling  is  that  the  tunneling  rates  are  highly  dependent  on  the  tunneling  probability, 
which  is  an  exponential  function  of  position  and  distance.  One  approximation  to  take 
is  to  assume  that  the  traps  at  each  oxide  node  are  most  likely  to  interact  with  their 
closest  interface  node.  This  would  require  a mapping  table  to  inform  the  simulator 
which  interface  nodes  should  be  linked  with  the  oxide  nodes  that  are  considered.  The 
components  of  the  Jacobian  matrix  can  then  be  modified  as  follows 


dF(nox)  _ _ _ dF(nox ) 

dnox  dnitf 


dF[nitf)  ' dF[ni,f) 

dnox  ' ' ' dnltf 


(4.22) 


4.3  Verification  with  Analytical  Expressions 

To  illustrate  the  oxide  trapping  noise  simulation  package  in  FLOODS,  a 0.4  gm 
n-channel  MOSFET  structure  with  a gate  oxide  thickness  of  140  A is  chosen  for 
testing.  The  noise  simulation  is  compared  with  the  analytical  expression  for  1 // 
noise  in  MOSFETs  to  determine  the  accuracy  of  this  tool.  In  FLOODS,  the  oxide 
defect  density  of  any  testing  device  structure  can  be  specified  in  the  form  of  Eq.  (4.18) 
proposed  by  Celik  and  Hsiang  [14].  Thus,  if  FLOODS  input  parameters  for  noise 
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calculation  are  properly  adjusted  to  include  only  the  physics  that  was  considered  in 
Celik  and  Hsiang’s  derivation,  comparisons  between  the  FLOODS  simulation  outcome 
with  their  theoretical  expectation  can  then  be  judged.  Under  the  assumption  of  an 
oxide  defect  density  Nt  ( E,x ) = Afoexp^  (£’  — £',•)]  and  with  the  square  barrier  tunneling 
model  to  characterize  the  transition  of  carriers  between  the  traps  at  the  interface  and 
in  the  oxide,  the  current  noise  spectral  density  at  the  drain  terminal  can  be  expressed 
as  [18] 


_ WjiWkTNo  f n0 SV 

S,d  21? fay  UJ  Vd(2Vd'sat  Vd> 
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{vd,saI-vdy 


vy 

d,sat 


(4.23) 


where  Vd  sat  = Vgs  — Vj,  y = £,kT,  and  hqs  is  the  surface  electron  concentration  at  the 
source.  Sid  of  Eq.  (4.23)  is  plotted  in  Figure  4-2  against  the  simulation  results,  and 
excellent  agreement  is  shown  in  the  figure  for  the  data  in  the  linear  region. 

In  the  past  the  1//  noise  was  explained  analytically  for  the  MOS  device  operating 
in  the  linear  mode  of  operation  only.  The  transfer  function  for  coupling  the  noise  from 
each  noise  source  was  derived  from  the  gradual  channel  approximation  (GCA)  and 
the  derivation  was  simplified  to  a one-dimensional  problem.  Also,  strong  inversion  is 
assumed.  Therefore,  the  1//  noise  formulation  breaks  down  in  saturation  and  in  sub- 
threshold operating  regions.  This  is  not  the  case  when  it  is  simulated  in  the  PDE-based 
simulator,  where  the  proper  transfer  function  can  be  computed  numerically  at  each  grid 
point,  in  either  the  sub-threshold,  linear,  or  saturation  modes  of  operation.  At  low  drain 
bias,  noise  simulations  are  performed  with  a frequency  sweep.  The  results  are  plotted 
in  Figure  4-3.  As  expected,  oxide  trapping  noise  is  a pure  1//  noise  with  a comer 
frequency  in  the  MHz  range.  Taking  a cross-section  of  the  distributed  current  noise 
spectral  densities  of  the  oxide  trapping  noise  at  a point  half  way  between  the  drain 
and  source  regions,  data  obtained  at  each  frequency  point,  a decade  apart,  is  plotted 
in  Figure  4-4.  It  is  shown  in  this  figure  that  the  first  line  from  the  left  corresponding 
to  the  lowest  frequency  (1  Hz)  peaks  around  22  A away  from  the  interface.  At  the 
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next  frequency  (10  Hz),  the  overall  magnitude  of  the  current  noise  spectral  density 
decreases,  and  at  the  same  time  the  peak  moves  closer  to  the  interface.  With  this 
figure,  one  can  identify  the  contributions  of  the  noise  produced  in  each  region  of  the 
oxide  at  any  specific  frequency,  thus  one  is  able  to  identify  the  critical  noise  region 
and  how  far  it  is  away  from  the  interface.  An  example  of  the  distributed  current  noise 
spectral  density  is  shown  in  Figure  4-5. 

To  observe  the  changes  in  oxide  trapping  noise  in  different  regions  of  device 
operation,  the  distributed  current  noise  spectral  density  in  the  oxide  is  plotted  under 
different  drain  bias  conditions.  As  stated  previously,  a square  barrier  tunneling  model 
with  a defect  density  distribution  independent  of  position  is  used  for  testing.  In  this 
model  the  tunneling  rates  and  the  characteristic  time  depend  only  on  distance.  This 
allows  the  peak  of  the  distributed  current  noise  spectral  density  to  be  located  at  the 
same  distance  away  from  the  interface  at  each  specific  frequency.  Once  the  critical 
distance  is  identified  for  the  chosen  frequency,  the  distributed  current  noise  spectral 
density  of  a cross-section  along  the  channel  under  different  drain  bias  is  plotted, 
as  shown  in  Figure  4-5  and  Figure  4-6.  Under  the  linear  mode  of  operation,  noise 
contributions  at  each  region  along  the  channel  increase  as  increases.  The  area 
under  the  curve  reaches  the  maximum  when  the  drain  bias  Vds  equals  Vd^at-  As  the 
drain  bias  is  increased  further,  the  peak  of  the  ridge  shifts  toward  the  source,  the  same 
way  the  channel  pinch  off  region  relocates.  Yet  the  area  under  the  curve,  which  is  a 
measure  for  the  noise  observed  at  the  contacts,  stays  the  same  when  compared  with  the 
curve  corresponding  to  Vds  = Vd,sat. 

No  significant  difference  to  noise  performance  was  found  by  switching  the 
tunneling  mechanism  from  the  square  barrier  tunneling  model  to  the  trapezoidal 
tunneling  model.  The  main  reason  for  such  an  insignificant  change  is  that  for  this 
device  with  a gate  oxide  thickness  of  140  A,  with  a gate  bias  of  1.0  V,  the  effect  of 
band  bending  at  the  locations  where  low  frequency  noise  is  significant  is  quite  small. 
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Thus  the  trapezoidal  barrier  seen  by  the  tunneling  carriers  looks  almost  like  a square 
barrier. 

As  mentioned  in  the  previous  chapter,  the  interface  g-r  noise  source  is  negligible 
when  compared  with  the  velocity  fluctuation  or  the  oxide  trapping  noise.  Thus  the 
interface  parameters,  such  as  the  interface  capture  coefficients  cn  and  cp  and  the 
interface  defect  density  Nritf,  cannot  be  extracted  from  noise  measurement  data.  In 
addition,  the  interface  trap  centers  serve  only  as  stepping  stones  for  the  carriers  in  the 
channel  to  reach  the  trap  centers  in  the  oxide,  and  the  carrier  transition  is  dominated 
by  the  carrier  tunneling  rates  gox  and  rox  defined  in  Equations  (4.13)  and  (4.14).  Thus, 
changing  the  values  of  the  interface  parameters  c„  and  cp  and  Njilf  does  not  alter 
the  simulated  oxide  trapping  noise  outcome,  as  long  as  those  parameters  ensure  that 
the  carrier  transition  rates  rn^t y,  gn  it y,  rpit y,  and  gpjt y are  always  higher  than  the 
tunneling  rates  gox  and  rox.  In  the  example  shown  in  this  section,  we  used  10-7  cm3/s 
and  1011  cm~2eV_1  for  the  interface  capture  coefficients  and  the  interface  defect 
density,  respectively.  The  parameter  which  determines  the  curvature  of  the  U-shaped 
trap  distribution  in  energy  is  set  to  E,  = 3.1,  as  determined  from  the  result  of  the 
conductance  method  [14,  18]. 

4.4  Noise  Spectroscopy 

Motorola  supplied  us  with  the  device  structure,  doping  profile,  and  measurement 
data  of  a 0.45  gm  graded  channel  MOS  (GCMOS)  transistor.  In  this  device,  a boron- 
doped  implant  is  present  near  the  source  region,  as  shown  in  Figure  4-7.  Both  the 
device  structure  (such  as  the  non-planner  oxide  thickness)  and  the  doping  profile  were 
calibrated  through  inverse  engineering  from  the  measured  I-V  and  C-V  data.  Once  the 
initial  device  structure  and  doping  profile  from  the  calibration  were  obtained,  all  oxide 
regions  of  the  device  were  re-gridded  with  fine  mesh  grids.  A prerequisite  for  accurate 
oxide  trapping  noise  simulation  is  to  have  dense  grid  points  in  the  oxide  meshes  near 
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the  interface,  hence  the  characteristics  of  the  noise  parameters  that  vary  strongly  with 
distance,  such  as  the  carrier  tunneling  rate,  can  be  accounted  for  properly. 

As  mentioned  previously,  the  interface  parameters  do  not  effect  the  oxide  trapping 
noise  outcome,  as  long  as  the  interface  carrier  transition  rates  are  faster  than  the  carrier 
tunneling  rates.  Thus,  for  Motorola’s  GCMOS  transistor  the  interface  parameters  for 
simulating  the  oxide  trapping  noise  are  set  to  cn  = cp  = 10-7  cm3/s  and  Njitf  = 1011 
cm  eV  . The  defect  density  in  the  oxide  is  modeled  according  to  Equation  (4.18), 
where  the  parameters  are  set  to  t,  = A = 3.1  and  r)  = 0.  For  the  carrier  tunneling  rate, 
the  attenuation  constant  a is  modeled  by  Equation  (4.21)  instead  of  (4.20)  for  better 
accuracy. 

Measured  noise  is  plotted  in  Figure  4-8.  The  data  was  taken  when  the  device 
was  biased  with  a constant  = 2.0  V and  sweeping  Vgs  from  0.4  to  0.9  V.  Since  the 
threshold  voltage  V,  extracted  from  both  the  measured  and  simulated  characteristic 
I-V  data  is  0.566V,  the  measurement  has  covered  both  the  subthreshold  and  inversion 
modes  of  operation  of  the  device.  To  analyze  the  excess  noise,  each  set  of  data  for  any 
specific  gate  bias  voltage  is  decomposed  into  two  parts:  the  background  pure  1//  noise 
component,  and  the  visible  g-r  Lorentzian  spectra.  We  speculate  that  the  Lorentzian 
spectra  visible  above  the  background  1//  noise  stem  from  discrete  defects  in  the  oxide. 
The  background  pure  1//  noise  for  each  gate  bias  is  represented  by  a solid  line  in 
Figure  4-8. 

To  analyze  the  background  1//  component,  simulations  were  carried  out  to  find 
the  constant  Nq  that  could  produce  the  noise  which  has  the  best  fit  to  the  measured 
data.  We  discovered  that  for  this  device,  the  critical  noise  region  is  not  in  the  gate 
oxide  above  the  center  or  near  the  drain  region  of  the  channel  like  the  example  shown 
in  the  previous  section,  but  is  in  the  smile  oxide  directly  above  the  graded  channel 
boron-doped  implant  near  the  source  region,  as  the  distributed  current  noise  spectral 
density  plotted  in  Figure  4-9  illustrates.  Even  though  the  peak  of  the  distributed  noise 
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moves  horizontally  (parallel  to  the  y axis)  with  the  changing  drain  bias  voltage  and 
moves  vertically  (parallel  to  the  x axis)  with  different  frequency,  the  peak  is  always 
confined  to  the  oxide  area  above  the  graded  channel  doping,  indicated  in  Figure  4-7. 
This  phenomenon  can  be  explained  by  a strong  scalar  Green  function  from  the  noise 
critical  region  to  the  drain  contact,  due  to  the  presence  of  the  graded  channel  doping. 

Studying  the  results  of  the  simulations  led  to  the  conclusion  that  the  parameter  No 
of  oxide  defect  modeling  in  Eq.  (4. 1 8)  cannot  be  limited  to  a constant  value.  Thus  we 
modeled  it  as  a function  of  the  position  y.  In  order  to  get  the  simulation  data  to  match 
the  measured  noise  data  for  a wide  range  of  gate  bias  values,  as  shown  in  Figure  4-10, 
the  parameter  Nq  in  the  oxide  noise  critical  region  needs  to  be  modeled  by  the  function 
shown  in  Figure  4-11.  This  result  demonstrates  that  the  defect  density  in  the  smile 
oxide  is  strongly  related  to  its  thickness.  As  derived  from  the  measured  noise  data,  the 
defect  density  is  higher  when  the  smile  oxide  is  thicker,  as  illustrated  in  Figure  4-12. 

Figure  4-10  shows  that  there  is  a discrepancy  between  the  simulated  and  mea- 
surement data  for  Vg.  = 0.4  V.  We  suspect  that  the  measured  noise  at  that  bias  is  not 
the  true  device  noise,  but  is  influenced  by  the  noise  introduced  from  the  measurement 
setup.  As  the  simulation  indicates,  the  critical  noise  region  for  the  gate  biases  Vgs  = 

0.4  V and  Vgs  = 0.5  V are  close  enough,  that  we  expect  that  the  same  discrete  defect 
centers  causing  the  Lorentzian  spectra  observed  in  the  noise  for  Vgs  = 0.5  V would  also 
show  up  in  the  noise  for  Vgs  = 0.4  V.  But  an  elevated  pure  1 If  noise  component  was 
observed  instead,  as  shown  in  Figure  4-8. 

To  observe  the  Lorentzian  spectra  in  the  measured  noise,  the  function  5/  x / 
vs.  / is  plotted  in  Figure  4-13.  From  this  figure  we  identify  three  major  Lorentzian 
spectra,  marked  by  a,  b,  and  c,  that  spread  among  the  five  sets  of  measured  data. 

To  approximate  these  spectra,  three  discrete  defect  centers,  a,  /3,  and  y,  are  placed 
inside  the  oxide  noise  critical  region,  shown  in  Figure  4-14,  to  produce  the  Lorentzian 
spectra  at  a,  b,  and  c in  the  spectral  window,  respectively.  By  adjusting  each  defect 
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center’s  location  and  magnitude,  we  obtain  a reasonably  good  agreement  between 
experiment  and  simulation,  as  illustrated  in  Figure  4-13.  A problem  with  the  discrete 
defect  distributions  is  that  they  may  spread  into  more  data  sets  than  was  intended.  For 
example,  the  discrete  defect  distribution  a,  which  intends  to  add  a Lorentzian  spectrum 
near  102  Hz  for  noise  measured  at  Vgs  = 0.8~0.9  V,  has  unintentionally  elevated  the 
simulated  noise  at  the  same  frequency  at  lower  gate  biases.  With  denser  grid  spacing 
and  more  iterative  steps  to  model  the  shape  of  the  discrete  defect  profiles,  one  can 
achieve  a better  match  between  the  simulated  and  measurement  data. 

4.5  Summary 

The  advantage  of  implementing  the  carrier-tunneling  mechanism  and  oxide 
trapping  noise  into  a PDE-based  simulator  is  that  one  is  able  to  observe  how  the 
l//-like  surface  noise  sources  are  distributed  in  the  oxide  and  at  the  interface.  Thus 
one  can  pinpoint  the  critical  noise  producing  region.  In  addition,  the  relationship 
between  device  surface  quality  and  the  l//-like  noise  can  be  simulated.  Comparing  the 
measured  and  simulated  noise  data,  oxide  and  interface  defect  density  can  be  estimated 
accurately.  This  information  can  be  used  to  predict  device  lifetime  and  estimate 
leakage  current  through  the  gate  oxide.  The  oxide  trapping  noise  simulation  can  also 
be  used  to  analyze  1/ /-like  noise  for  MOS  devices  operating  in  the  saturation  mode. 
This  cannot  be  done  analytically,  due  to  the  2-dimensional  complexity  of  the  channel 
layout  and  carrier  flow  patterns. 


Figure  4-1.  “Tunneling  and  Capture”  and  “Capture  and  Tunneling”  models 
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Figure  4-2.  Simulation  results  are  plotted  and  compared  with  the  analytical  prediction. 
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Figure  4-3.  Current  noise  spectral  density  vs.  frequency 
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Figure  4-4.  Distributed  current  noise  spectral  density  for  different  frequencies 
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Figure  4-5.  Distributed  current  noise  spectral  density  along  the  channel  for  Vgs  = 1.0  V 
and  Vds  = 0.45  V (/  = 1 Hz) 
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Figure  4-6.  Distributed  current  noise  spectral  density  along  the  channel  for  different 
bias  conditions 
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Figure  4-7.  Device  structure  of  a 0.45  fim  graded  channel  MOS  transistor 


Noise  Current  Spectrum  (A  /Hz-pm) 
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Figure  4-8.  Measured  current  noise  spectral  density  as  a function  of  Vgs  and  frequency. 

The  drain  bias  Vjs  = 2.0  V.  The  solid  lines  represent  the  simulated  back- 
ground pure  1 If  noise  components. 
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Figure  4-9.  Distributed  current  noise  spectral  density  for  / = 25  Hz.  The  bias  values 
of  the  device  under  simulation  are  = 2.0  V and  Vgs  = 0.8  V.  The  oxide 
region  is  shown  with  x = 0 to  be  the  midway  between  the  metallurgical 
source  and  drain  junctions. 


Noise  Current  Spectrum  (A  /Hz-jim) 
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Figure  4-10.  Comparison  between  the  simulated  and  measured  data  of  the  pure  1 // 
noise  component  (V^  = 2.0  V,  / = 25  Hz). 
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Figure  4-11.  Profile  of  the  defect  density  parameter  Nq  as  function  of  position  parallel 
to  the  channel.  Note  that  only  the  parameter  in  the  oxide  noise  critical 
region  indicated  in  Figure  4-7  is  shown  here. 


88 


Figure  4-12.  Defect  density  profile  in  the  smile  oxide 
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Frequency  (Hz) 


Figure  4-13.  Product  of  current  noise  spectral  density  and  frequency  to  identify 
Lorentzian  spectra  in  the  measured  and  simulation  data 
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Figure  4-14.  Discrete  defect  distributions  responsible  for  the  Lorentzian  spectra  ob- 
served in  the  experiment.  Note  that  only  the  distributions  in  the  oxide 
noise  critical  region  indicated  in  Figure  4-7  are  shown  here. 


CHAPTER  5 

MAXIMUM  ALLOWABLE  BULK  DEFECT  DENSITY  FOR 
GENERATION-RECOMBINATION  NOISE  FREE  OPERATION  OF  n+-n-n+ 

STRUCTURES 

5.1  Introduction 

Many  defects  that  exist  in  semiconductor  devices  are  often  unintentional.  Elements 
such  as  Carbon,  Oxygen,  and  various  metals  may  be  introduced  to  the  wafer  during 
processing  steps,  and  then  become  carrier  trap  centers  that  affect  device  performance 
and  produce  generation-recombination  (g-r)  noise.  The  g-r  noise  is  due  to  the  ran- 
dom fluctuations  in  the  trapping  and  de-trapping  rates  in  the  defect  trap  centers  in  the 
forbidden  band  gap.  This  in  turn  causes  the  current  or  the  voltage,  depending  on  the 
condition  at  the  contact  terminal,  to  fluctuate  as  well.  The  generation  and  recombi- 
nation transition  has  generally  been  characterized  by  the  Shockley-Read-Hall  model. 

In  this  model,  each  defect  species  is  specified  by  its  defect  density,  its  carrier  capture 
coefficients,  and  its  energy  position  relative  to  either  the  conductance  or  valence  band. 
As  has  been  shown  in  the  past,  g-r  noise  offers  a unique  opportunity  to  experimentally 
determine  these  trap  parameters  on  as-made  devices  under  actual  bias  conditions. 

An  advantage  of  the  g-r  noise  trap  spectroscopy  over  other  defect  spectroscopy  tools 
is  that  only  defects  that  affect  device  performance,  via  their  noise  contributions,  are 
probed  [26].  The  focus  of  this  paper  however  is  to  study  the  maximum  allowable 
bulk  defect  density  to  guarantee  g-r  noise  free  operation  under  all  bias  conditions  of 
specifically  two  terminal  drift  dominated  devices  which  are  known  to  be  sensitive  to 
bulk  defects.  Having  this  information  may  guide  the  decision  on  the  need  for  incor- 
porating defect  reduction  steps  in  the  device  fabrication  process.  Both  the  local  g-r 
noise  strength  and  the  transfer  Green’s  function,  which  couples  the  carrier  fluctuations 
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in  each  differential  volume  of  the  device  to  a specific  contact,  change  with  bias.  The 
simple  analytical  model  of  g-r  noise  indicates  that  the  current  noise  spectral  density 
of  a bulk  resistor  in  the  ohmic  regime  is  proportional  to  72  or  E2,  where  / is  the  cur- 
rent and  E is  the  electric  field  in  the  device.  As  a consequence,  increasing  the  bias 
across  a device  will  nearly  always  lift  a g-r  noise  component  out  of  the  thermal  noise 
background.  In  practice  however,  carrier  heating  and  space  charge  injection  from  the 
contacts  will  lead  to  field  dependent  carrier  mobilities  and  position  dependent  carrier 
densities  resulting  in  a decrease  of  the  g-r  noise  signature  of  a defect  center  at  high 
bias  to  values  possibly  below  the  velocity  noise  level.  It  is  difficult  to  unravel  the  full 
physical  noise  picture  analytically.  We  therefore  expanded  the  capabilities  of  a semi- 
conductor device  simulator  with  g-r  and  velocity  fluctuation  noise  models  to  account 
for  these  effects. 


The  g-r  process  can  be  described  by  the  Shockley-Read  Hall  model,  as  shown  in 
Fig.  1.  For  illustration  purposes,  assume  there  is  one  defect  trap  level  in  the  bandgap, 
and  the  traps  are  acceptor-like.  In  addition  to  the  three  basic  Shockley  equations 


5.2  Implementation 


Fv/  = ~^~e  ip~n  + ND~NA  ~nt ] =0 
dn  1 


(5.1) 


(5.2) 


(5.3) 


a trapped  electron  continuity  equation  is  added  to  describe  the  carrier  density  fluctua- 
tions at  the  trap  level  Ej 


(5.4) 


where  y„,  yp,  and  y are  the  Langevin  noise  terms  describing  random  transition  rate 
fluctuations.  The  transition  rates  can  be  expressed  as 
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rn  — cnn{NT  — nt) 


(5.5) 


8n  — cnn\tit 


(5.6) 


rp  = cppnt 


(5.7) 


gP  = cPP\  (Nr-n,) 


(5.8) 


where  Nj  represents  the  defect  density  and  n\  and  p\  are  the  Shockley  parameters 
defined  as 


with  the  assumption  that  the  trap  level  is  at  least  a few  kT  from  the  conductance  and 
valence  band  edges  and  other  symbols  have  their  usual  meaning.  In  these  equations  the 
defect  species  are  specified  by  three  parameters:  trap  energy  level  (Et),  trap  density 
(Nr),  and  capture  coefficients  (cn,  cp). 

Physics  based  noise  models  are  implemented  into  the  device  simulator  using  the 
impedance  field  method  (IFM)  [4].  Two  elements  are  required  to  calculate  the  noise 
at  the  contact  terminals.  The  first  is  the  magnitude  of  the  microscopic  noise  sources 
Kya,yp  (where  a,/3  = n.p,nt)  at  each  grid  point,  represented  by  the  added  Langevin 
noise  terms  in  the  partial  differential  equations  (PDE)  in  (5.2)-(5.4).  For  the  g-r  noise 
sources,  the  microscopic  noise  sources  Kya ^ are  derived  from  the  Poisson  process  and 
are  given  by  [11] 


(5.10) 


(5.9) 


Kyn,yn  — 2gn  + 2rn  = — £y„,y, 
Kyp,yp  — 2gp  + 2 rp  = —Kypiyt 
Kyt,y,  = Kyn,yn+KyP  ,yP 


(5.11) 
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The  second  element  is  the  transfer  function  that  couples  the  perturbation  of  the  PDEs 
from  each  grid  point  to  the  device  terminals.  The  scalar  Green’s  transfer  function  G for 
each  carrier  type  at  each  node  is  obtained  by  perturbing  the  PDE  associated  with  the 
earner  at  the  node,  and  then  observing  the  changes  at  the  contact  terminals.  The  g-r 
noise  current  spectral  density  at  an  external  contact  terminal  is  calculated  through  the 
following  expression  [4] 


where  Ntrap  is  the  number  of  trap  levels  in  the  device. 

The  g-r  noise  mechanism  described  above  was  implemented  into  FLOODS 
(Florida  Object-Oriented  Device  Simulator).  FLOODS  is  a PDE-based,  generalized  box 
scheme  device  simulator.  The  implementation  of  velocity  fluctuation  noise  is  described 
in  an  earlier  paper  [6]. 


Extensive  research  on  bulk  g-r  noise  for  either  single  or  multiple  defect  trap 
systems  [11,  12]  has  shown  that  the  noise  spectrum  of  the  number  of  trapped  carriers 
Nt  in  a small  differential  volume  Av  equals 


(5.12) 


5.3  Simulation 


(5.13) 


where  co  = 2nf.  The  spectrum  of  yt  is  equal  to 


Sy  — lim  2iy;)f  — 2(gn  + rn+gp  + rp) 

T — >°° 


(5.14) 


where  the  units  of  yt,  and  the  generation  and  recombination  rates  have  been  changed 
into  number  of  carriers  per  second  for  the  remainder  of  this  paper.  Combining 
(5.13,5.14)  results  in 


(5.15) 
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The  value  of  the  low  frequency  plateau  is 

San,,Av  (f  0)  — 2 (gn  + rn+  gp  + r p ) T2 


(5.16) 


To  simulate  bulk  defect  noise  properties,  a 2-D  n+-n-n+  bulk  resistor  was  created.  The 
n-  and  n+-region  are  doped  with  Nj  = 1015  and  Nj+  = 1018  cm-3  donors,  respectively. 
The  length  of  the  n-region  L is  2 /im,  and  the  n+-regions  are  0.5  pm  long  in  length  on 
each  side  of  the  device.  The  cross-sectional  area  A is  1 pm2.  A defect  center  Ep  at  Ep 
= ( Ec  - 0.265  eV)  is  chosen  for  the  simulation.  This  choice  was  guided  by  the  fact  that 
trap  levels  close  to  Ep  produce  maximum  noise  levels  [6]  and  thus  represent  a worst 
case  scenario.  Electron  and  hole  capture  coefficients,  c„  and  cp  and  defect  density  Np 
are  set  to  10~8  cm3/s  (typical  values)  and  , initially,  1012  cm-3,  respectively.  Since 
Ep  is  much  closer  to  Ec  than  Ey,  one  can  safely  assume  gn,rn»  gp,rp  and  g„  « rn. 
Eq.  (5.15)  above  then  simplifies  to 


With 


San,,Av  — 


4rnz? 

1 + £02t,2 


4g„T2 

1 + CD2  X} 


(5.17) 


San,  - San 


(5.18) 


where  AN  is  the  number  of  carriers  in  the  conduction  band  in  the  differential  volume 
and  integrating  (5.17)  over  device  volume,  the  current  spectral  density  of  the  resistor  in 
the  ohmic  regime  becomes 


Sl  (/)  = ]y2 ' Sw  (/) 


4/2g„T2  1 

(LAn)1  1 + «2t2 


(5.19) 


where  we  used 


Sn  _ Sy  _ Si 
N2  ~ V2  ~ 1 2 


(5.20) 


The  symbols  V and  I are  the  dc  bias  voltage  and  current,  respectively,  and  N represents 
the  total  number  of  carriers  in  the  device.  Since  cn(n  + n\)  » cp(p  + p\),  xt  can  be 
simplified  to  xt  « [cn  (w  + ni)]_1.  The  analytical  expression  of  the  current  noise  spectral 
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density  of  this  structure  at  low  frequencies  becomes  then 


(5.21) 


The  g-r  noise  of  the  n+-n-n+  bulk  resistor  with  the  defect  trap  center  described 


above  is  simulated  and  plotted  as  a function  of  frequency  in  Fig.  5-2  and  compared 
with  the  analytical  g-r  noise  expression  (5.19).  The  simulated  noise  is  in  excellent 
agreement  with  the  theory.  The  simulated  bulk  g-r  noise  rolls  off  at  the  calculated 
comer  frequency  fc  = \/2nzt  « 3 MHz,  and  has  a slope  of  l//2  beyond  fc. 

G-r  noise  simulations  as  a function  of  bias  are  shown  in  Fig.  5-3.  Under  low 
bias  conditions,  the  g-r  noise  is  below  the  velocity  fluctuation  noise.  The  simulation 
is  in  good  agreement  with  the  analytical  expressions  (5.19)  and  (5.21)  under  low 
bias  conditions.  When  the  bias  increases,  the  magnitude  of  the  simulated  g-r  noise 
emerges  from  the  velocity  fluctuation  noise,  reaches  a peak  and  then  starts  to  decrease 
due  to  predominantly  space-charge  injection  from  the  n+  drain  region  of  the  device. 
The  space  charge  injection  as  a function  of  bias  is  illustrated  in  Fig.  5-4.  Due  to  the 
injection,  injected  electrons  fill  the  traps  with  a doping  density  of  1012  cm-3  as  shown 
in  Fig.  5-5  and  subsequently  the  electron  Fermi  level  in  the  n-region  moves  away  from 
the  trap  energy  level.  As  a result  the  g-r  mechanism  at  that  level  becomes  less  effective 
due  to  a decrease  in  the  recombination  rate  and,  via  detailed  balance,  generation  rate. 
The  result  is  then  a decrease  of  the  local  g-r  noise  strength. 

The  scalar  Green’s  function  can  be  derived  from  (5.20)  analytically,  as 


(5.22) 


Under  low  bias  condition  where  N is  independent  of  bias,  the  Green’s  function  is 
proportional  to  I.  Under  high  bias  conditions,  N increases  with  bias  as  the  carrier 
injection  into  the  n-region  becomes  significant.  As  a result  the  Green’s  function 
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becomes  less  bias  dependent,  as  shown  in  Fig.  5-6.  The  same  trend  can  be  seen  in  the 
distributed  current  spectral  density,  shown  in  Fig.  5-7. 

5.4  Discussion 

Fig.  5-3  shows  that  for  the  trap  example  selected  the  g-r  noise  dominates  the 
velocity  fluctuation  noise  over  a wide  bias  range,  hence  the  chosen  trap  density  of 
1012  cm-3  clearly  exceeds  the  maximum  allowable  defect  density  NTmax  for  g-r  noise 
free  operation.  The  value  of  Njmax  for  this  particular  trap  follows  from  equating  the 
g-r  noise  magnitude  with  the  magnitude  of  the  velocity  fluctuation  noise  component. 
Under  low  bias  conditions  in  the  ohmic  regime,  the  quantity  of  Njmax  can  be  derived 
analytically  by  setting  (5.21)  equal  to  S/av  = 4kT/R.  Thus 

kTL 2 cn(n-\-n\)i 


NTn 


q/dV 2 


n\ 


(5.23) 


with  the  assumptions  that  nt  ~ c„nNtzt  and  zt  ~ [cn  (n  + «i)]_1.  Simulations  were 
carried  out  by  varying  the  value  of  Nt  to  obtain  a g-r  noise  level  equal  to  the  velocity 
fluctuation  noise  level  as  a function  of  trap  energy  level  for  a bias  of  0.1  V at  room 
temperature  with  the  c„  and  cp  values  previously  specified.  The  simulation  results  are 
presented  in  Fig.  5-8  and  agree  with  the  low  bias  predictions  of  (5.25).  The  figure 
shows  that  defect  levels  close  to  the  Fermi  level  at  Eq  — Ef  = 0.265  eV  require  stricter 
constraints  than  defect  levels  further  away  from  Ep , corroborating  the  fact  that  traps  at 
the  Fermi  level  produce  g-r  noise  more  effectively.  Fig.  5-8  also  shows  that  at  higher 
bias  lower  defect  densities  are  required  for  g-r  noise  free  operation  primarily  due  to  an 
increase  in  magnitude  of  the  Green’s  function  with  bias.  To  estimate  NTmax  at  high  bias 
for  a trap  at  the  Fermi  level,  the  worst  case  scenario,  we  use  n\  = n and  n,  = 0.5  Nt. 
Upon  substituting  these  assumptions  in  (5.17)  and  integrating  over  device  volume,  the 
low  frequency  plateau  value  of  the  number  spectral  density  becomes 


NtAL 
2 cnti 


San{Q) 


(5.24) 
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The  value  of  San  (0)  via  (5.22)  needs  to  be  compared  with  the  magnitude  of  the 
velocity  fluctuation  noise  S/av  component  as  read  from  the  simulated  data  presented 
in  Fig.  5-3.  Note  that  I,  N,  and  S/av  are  affected  by  carrier  heating  and  space  charge 
injection.  The  value  of  I follows  from  the  simulated  I-V  characteristic  presented  in 
Fig.  5-9,  whereas  N follows  from 


CV 

N=  — +NdAL 

q 


(5.25) 


where  C = \.1CQ  with  C0  representing  the  geometrical  capacitance  of  the  device. 

Eq.  (5.25)  is  based  on  early  work  by  Lampert  and  Mark  [27]  and  was  later  confirmed 
by  Bosman  et  al.  [28],  Values  of  N thus  calculated  agree  well  with  simulated  values 
presented  in  Fig.  5-4.  The  value  of  Njmax  can  now  be  calculated  from 


NTmax  = 2 


cnN3SIAv(  0) 


(ALI)2 

resulting  in  Njmax  = 1.4-1011  cm-3  at  2 V and  2-1011  cm-3  at  4 V,  respectively,  in 
very  good  agreement  with  the  simulated  data  presented  in  Fig.  5-8. 

5.5  Conclusions 

A detailed  comparison  of  simulated  and  analytical  noise  predictions  shows  that 
the  maximum  allowable  defect  density  for  g-r  noise  free  operation  of  two  terminal 
bulk  devices  can  accurately  be  estimated  even  in  the  presence  of  hot  carrier  effects 
and  space  charge  injection.  In  addition  the  noise  simulations  reveal  that  the  g-r 
noise  plateau  value  of  a specific  trap  reaches  a maximum  as  a function  of  bias  and 
subsequently  decreases  with  increasing  bias  due  to  space  charge  injection  and  carrier 
heating  to  levels  below  the  velocity  fluctuation  noise. 


(5.26) 
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Figure  5-1.  Carrier  transitions  between  the  continuum  states  and  localized  defect  states 


Noise  Current  Spectrum  (A2/Hz) 


100 


Figure  5-2. 


Comparison  between  the  simulated  and  analytical  current  noise  spectral 
density  as  a function  of  frequency.  The  bias  across  the  device  is  set  to 


0.1  V 


Noise  Current  Spectrum  (A  /Hz) 
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Figure  5-3.  Comparison  between  the  simulated  and  analytical  g-r  current  noise  spectral 
density  (/  = 1 Hz)  as  a function  of  the  bias  voltage.  The  dashed  line  rep- 
resents the  analytical  expression  (5.21).  The  solid  line  is  the  simulated  g-r 
noise.  The  dotted  line  is  the  simulated  velocity  fluctuation  noise. 


Electron  Concentration  (cm 
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Figure  5-4.  Electron  concentration  n in  the  n-region  as  a function  of  position  under 
different  bias  conditions.  Under  the  low  bias  condition  (V^  = 0.1  V), 
nmNd.  As  the  bias  increases  (VbUu  > 1.0  V),  n > Nj. 
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Figure  5-5.  Trap  electron  concentration  nt  as  a function  of  position  under  different  bias 
conditions 


Real  part  of  scalar  Green’s  Function  (As) 
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Figure  5-6.  Scalar  Green’s  function  under  different  bias  conditions 


Distributed  Noise  Current  Spectrum  (A2/Hz-cm3) 
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Figure  5-7.  Distributed  current  noise  spectral  density  under  different  bias  conditions 
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Figure  5-8.  Defect  density  Njmax  as  a function  of  defect  energy  under  different  bias 
conditions 


Current  (A) 
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Figure  5-9.  I-V  characteristic 


CHAPTER  6 

CONCLUSION  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH 

The  velocity  fluctuation  noise,  the  bulk  and  interface  generation-recombination 
(g-r)  noise,  and  the  oxide  trapping  noise  mechanisms  were  implemented  into  FLOODS. 
Each  has  been  verified  with  its  known  analytical  expressions  to  validate  its  accuracy. 

By  matching  the  simulation  with  the  measured  noise  data  of  the  device,  this  simulation 
tool  can  become  useful  for  probing  the  defect  density  distributions  in  both  the  silicon 
bulk  and  oxide  regions. 

One  suggestion  for  future  research  is  to  utilize  the  noise  simulation  tool  to  study 
the  noise  in  advance  VLSI  devices,  such  as  silicon-on-insulator  (SOI)  MOS  transistors. 
The  carrier  tunneling  mechanism  that  produces  the  1//  noise  in  MOS  devices,  as 
the  author  has  demonstrated,  also  exists  in  SOI  devices.  It  is  possible  that  the  oxide 
trapping  noise  is  more  critical  for  SOI  devices  than  for  MOS  devices,  since  they  have 
more  silicon-oxide  interface  area.  In  addition,  since  the  carrier  tunneling  is  not  limited 
by  the  device  size,  the  oxide  trapping  noise  simulation  can  also  be  used  to  analyze  the 
low-frequency  noise  in  short-channel  (L  < 0.4  gm)  MOS  transistors. 

As  both  the  length  and  the  width  of  the  MOS  devices  become  smaller,  single 
defects  in  the  oxide  may  become  important.  Instead  of  getting  a l//-like  noise 
spectrum,  the  simulated  outcome  would  consist  of  a limited  numbers  of  Lorentzian 
spectra,  and  show  Random  Telegraph  Signal  (RTS)  noise  in  the  time  domain.  Since  the 
noise  simulation  tool  is  capable  of  accepting  discrete  defect  density  distributions,  it  can 
be  used  to  study  and  simulate  the  RTS  noise. 

Another  suggestion  is  to  combine  the  oxide  trapping  noise  with  the  harmonic 
balance  simulation,  to  simulate  the  1 If  noise  sidebands  observed  in  the  mixers  and 
voltage  controlled  oscillator  (VCO)  circuits. 
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Finally,  another  interesting  topic  would  be  to  use  the  oxide  trapping  noise  simu- 
lator to  investigate  the  mobility  fluctuation  noise.  As  the  carriers  tunnel  between  the 
traps  in  the  oxide  and  at  the  interface,  the  perpendicular  field  at  the  channel  would  fluc- 
tuate (due  to  the  charging  and  discharging  of  the  trap  centers  as  the  carriers  tunneling 
in  and  out).  This  indirectly  modulates  both  the  carrier  mobility  and  current  density  in 
the  channel.  Preliminary  simulation  results  have  shown  that  this  phenomenon  can  be 
captured  by  the  noise  simulation  tool.  More  research  is  currently  underway  to  study 
this  phenomenon. 
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